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ABSTRACT: The Sekine sensitivities to pressure and temperature of three t hal a 
- emitters (helical long-path, labyrinth and vortex-type emitters) were determined in the 4 
_ laboratory. The measured discharge sensitivity to pressure was used to calculate the 
- maximal permissible pressure variation in a field for the various emitters. The relative - 
_ discharge distribution along an experimental pipe with helical long-path emitters was 
_ calculated on the basis of the measured temperature distribution and the computed — Ty 
distribution. The resulting relative discharge distribution was curvilinear, in 

_ good agreement with the measured values. Under field conditions where water is - 
heated as it flows in pipes exposed to solar radiation, emitter discharge along the 
second half of the pipe is determined by its specific discharge sensitivity to | 


and Temperature,” Journal of the Irrigation and shite Division, ASCE, Vol. 107, a 
No. IR1, Proc. Paper 16116, March, 1981, 


Zur, and Tal, S., “Emitter Discharge Sensitivity Pressure 


16115 ‘HISTORY OF FEDERAL WATER RIGHTS 


KEY WORDS: Federal control; « Federal-s state. 
_ History; Water allocation (policy); Water law; Water te cooperation; | “7 


Water resources development; Water rights 

ABSTRACT: The historical development of Federal reserved water rights, of interest 
_ to all who are involved in water-resource planning or use, is presented. The Federal 
_ reserved rights are those that are reserved in fact or by implication in Federal treaties, 


| reservations, acts, or actions. Federal actions include flood-control projects and _ 


- navigation improvement. The courts of the United States have developed and 
_ promulgated a system of Federal reserved water rights that can be traced from about 
_ 1850 to existing suits on water resources. A determination of Federal water rights is 


necessary for state water-resource planning and allocation. sag 


_ REFERENCE: Bird, John W., “Origin and Growth of Federal Reserve Water Rights,” 


Journal of the Irrigation and Drainage Division, ASCE, Vol. 107, No. IR1, Proc. 


Paper 16115, March, 1981, pp. 11-24 
pp. 11-240 
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| 16101 AQUIFIER ESTIMATION AND KALMAN 


KEY WORDS: Aquifers; Filters; Groundwater sources; Parameters; Pump L 


H 
> oe ABSTRACT: An iterated extended Kalman filter (IEKF) has been used for the 


q 


ii estimation of aquifer parameters in the presence of uncertainties in the system and 

- noise in the measurements. The technique is sequential and each additional observation | 7 

the error covariance of ‘parameter estimates. The method is 
computationally efficient and gives the confidence limits of the parameter estimates. 

_ Only a few observations are sufficient to provide a reasonable estimate of the —~ il 
parameters. Actual aquifer test data for confined nonleaky and leaky aquifers —/ 

_ been analyzed and the results compared with those obtained using the known 

| The IEKF method gives less residual square error and eliminates the 
subjectivity in the conventional curve matching d 


*. REFERENCE: Chander, Subhash, Kapoor, Prakash N., and Goyal, Sushil K., “Aquifer 
Bi Parameter Estimation Using Kalman Filters,” Journal of the Irrigation and Drainage _ 
Division, ASCE, Vol. 107, ‘No. IR1, Proc. bees 16101, March, 1981, pp. 25-33 
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KEY WORDS: Groundwater chemistry; Gremiuater quality; Irrigation; — 
Irrigation effects; Leaching; Precipitation (chemistry); Salinity; Travel _ = Aele, 
"Increased irrigation efficiency usually reduces the degradation of 7 
 secaiadiennan’ resources. The effect of differences of travel times in the unsaturated zone _ 
between low and high leaching becomes more important as the depth to the 
groundwater table increases. Waters approaching gypsum saturation have a much . 7 
_ greater potential for ground-water degradation at high leaching as compared to low | 
leaching irrigation management. This is due to higher salt concentrations and, we 
|! consequently, increased epee Precipitation in the unsaturated zone in the case of | a 


1 REFERENCE: Suarez, Donald L., and van Genuchten, Martin Th., “Leaching and 
| Water-Type Effects on Ground-Water Quality,” Journal of the Irrigation and Drainage 


ao Division, ASCE, Vol. 107, No. IR1, Proc. Paper 16119, March, 1981, pp. 35-520 


KEY WORDS: Crop production; Crop yield; Economic factors; Irrigation; Bt: 
Models; River basin planning; River basins; Salinity; Seasonal variations; "7 Co 


‘method is presented that combines crop ‘model with 

a hydrologic-salinity simulation model to estimate irrigation management effects of — q 
sequential water use in a river basin. Crop yields are assumed to be a function of _ 
seasonal evapotranspiration and soil solution salinities. The combined model was _ 
applied to the Hadejia River Basin in northern Nigeria. The relative efficiency of water 
& with respect to crop production, is presented for selected management alternatives. 

REFERENCE: Fapohunda, Henry O., and Hill, Robert W., “River Basin Hydro-— 
Salinity-Economic Modeling,” Journal of the Irrigation and Drainage Division, ASCE, 
+ 107, No. IR1, Proc. Paper 16145, March, 1981, pp. 53-69 


16134 INFILTRATION. EQUATIONS FOR SURFACE EFFECTS 
KEY WORDS: Computer models; Equations; Hydrology; Infiltration 
capacity; Infiltration rate; Infiltration (soils); Mulching; Runoff; Seating: 
ABSTRACT: A series of equations were predicting infiltration into soils 
with modified surfaces. The Green-Ampt approach was used in. their development. 
_ They have the same functional form as the Green-Ampt and Mein-Larson two-stage — 
_ infiltration equations to which they reduce for the case of a uniform soil profile. The — 
- equations were applied to coarse over-fine and fine over-coarse soil stratifications and 
_ showed good agreement with both observed infiltration rates and infiltration rates 
_ predicted by numerical solution of Richards’ equation. The equations were also applied 
to an Ida silt loam soil to illustrate the magnitude of the predicted effect of surface _ 
_ Sealing on infiltration. Analysis was performed for two different initial moisture — 
contents, three different rainfall intensity ratios, and four different conductivity ratios. a 
q _ REFERENCE: Moore, lan D., “Infiltration Equations Modified for Surface Effects,” 
Journal! of the Irrigation and Drainage ASCE, Vol. 107, No. Proc. 
| Paper 16134, March, 1981, pp. 71-86 
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16139 HEAD LOSS OVER LONG: THROATED FLUMES 


Downstream; (fluid flow); Flumes; Laboratory 
“ABSTRACT: Long- flumes, all with th ame trapezoidal throat, 
laboratory tested. The flume throat was first fitted with two entrance transitions with 
plane bottoms converged respectively 1:2, and 1:3, as measured relative to the flume io 
center line. The 1:2 transition has no significant influence on the flume calibration. The 
‘throat was then tested with six downstream transitions of which the expansion ratio ’ 
ranged from 1:0—1:10. Two flume series were tested with these transitions. The . 
modular limits and related required energy head losses were analyzed. An equation is 
presented to estimate the modular limit for any long-throated flume placed in an 
arbitrary canal; a procedure to calculate this modular limit is also presented. 


_ REFERENCE: Bos, Marinus G., and Reinink, Yvonne, ‘ “Required Head Loss over 


Long- -Throated Flumes,” Journal of the Irrigation and ASCE, 

107, IRI, Proc. Paper 16139, March, 1981, PP. 87-102 
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; In accordance with the October, 1970 action of the ASCE Board of Direction, which stated 
t that all publications of the Society should list all measurements in both U.S. Customary and 
; ‘SI (International System) units, the following list contains conversion factors to enable readers 
0 compute the SI unit values of measurements. A complete guide to the SI system and its | 
e has been published by the American Society for Testing and Materials. Copies of this 
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publication (ASTM E-380) can be purchased from ASCE at a price of $3.00 each; orders must 
All authors of Journal papers are being asked to prepare their papers in this dual-unit format. 


“miles (miles) 


square inches (sq in.) 

‘square feet (sq ft) 7 ~ square meters (m i 

‘Square yards (sq yd) square meters (m*) 

square miles (sq mils) square kilometers (km?) 
acres (acre) hectares (ha) 
cubic inches (cu in. cubic millimeters (mm° 
cubic feet (cu ft) cubic meters (m 


= cubic meters (m 
tons (ton) mass kilograms (ke) 


kilogram force (kgf) newtons (N) 
pounds per square foot (psf) — pascals (Pa) “a 
U. S. gallons i 1 
The 


| 
| __kilometers (km) inl 
645 a 7 
0.093 
0.836 4 
90-0280 
= 
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Emrrrer D DISCHARGE SENSITIVITY TO PRESSURE 


as 


nied! an irrigation line and over the field. In the design of an irrigation network, 
emitter discharge is traditionally assumed to be a function of pressure only. — 
ad Discharge uniformity is assured by limiting line pressure variation over the 
field. The upper and lower limits of line pressure variation would depend on f 
the discharge sensitivity to pressure of the specific emitter. 
e _ ‘The temperature of water flowing in irrigation | pipes exposed to solar radiation 
be is expected to increase with distance. Parchomchuk (3) measured a 16° C (61° F) ; 
ers in the temperature of the water as it passed through a 37-m (121-ft) long” 
_ polyethylene pipe exposed to the sun. He also measured a 6° C (43° F) rise i 
4 in water temperature in a similar pipe buried 15 cm (6 in.) below the ssithnn: 
> Gilad et al. @ measured the water ‘temperature distribution along a 70-m (229-ft) _ 
a a 12°C (54° F) rise in the a 
3 3 water cemperature at the end of the line. The major part of the at as 
fise was measured along the last 15 m (49 ft) of the pipe. 
_ The sensitivity of emitter discharge to temperature was oy 
“and Karmeli (2). They suggest that the discharge sensitivity to temperature _ 
_ would be be minimal in in emitters where the flow regime is turbulent, e.g., a 
jortex- , and labyrinth- type emitters. For helical long- path emitters, Keller 
and Karmeli (2) assume that since the flow regime is laminar, the discharge — 
as Be senemavity to temperature will be equal to the sensitivity of viscosity to ane 


‘The objective of the present work was to measure the discharge sensitivity 
to pocesure and temperature for three types of emitters, and to investigate the | 
"Sr. Lect., 


Dept. of Agricultural Engrg., Technion, Haifa, Israel. ts” 
?Grad. Research Asst., Dept. of Agricultural Engrg., Technion, Haifa, Israel. cent 
Note.—Discussion open until August 1, 1981. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, 
‘ fi ASCE. Manuscript was submitted for review for possible publication on June 26, 1980. 
This paper is part of the Journal of the Irrigation and Drainage Division, Proceedings x 
of the American Society of Civil ©ASCE, No. March, 1981. 
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relative importance of pressure and temperature distribution on he discharge 


AGI 
Three types of emitters, by mechanism, 
P* were used in the present study. The helical long-path principle was represented 


; c by four ‘‘Netafim’’ emitters with nominal discharges of 1 L/h, 2 L/h, 4 L/h, 


_ and 8 L/h. The labyrinth principle was represented by the “‘Lego 4’’ emitter 
™ with a nominal discharge of 4 L/h, and two ‘‘Kochav’’ emitters with nominal 
discharges of L/h and 4 L/h. 1 The vortex Principle was by the 


4, Kochay, and Tirosh = 5 are all commercial names of emitters manufactured 
_ The pressure and temperature sensitivity | of emitter discharge was determined - 
for each of the aforementioned emitters. These tests were conducted using — 

160-cm long plastic pipes on which five emitters of the same kine 

at 30-cm intervals. During each test, the discharge of each emitter was measured 
three times. The pressure-discharge relationships of each kind of emitter were 
determined by measuring the discharge at six line pressures of 50 kPa, 100 
kPa, 150 kPa, 200 kPa, 250 kPa, and 300 kPa. The sensitivity of emitter discharge — 

to temperature was measured at seven water temperatures in the range of 

7° C-55° C (44° F- -131° F), and | at three line pressures of 100 kPa, 200 kPa, 
and 300 kPa. Water temperatures above ambient were obtained by mixing tap B 

water with hot water from a commercial water heater in various proportions. 

Water temperatures below ambient were obtained by passing tap water through 

a copper pipe immersed in a container filled with ice water. Water temperatures 


were € measured at the outlets of the first, ‘third, and fifth emitter ¢ on the | annul 


& 


Discharge to —The pressure- ure-discarg of 


which Q = = the the pressure; and a b= 
‘The values of a and b for each emitter were determined by a linear cosa. 
analysis of a logarithmic plot of the measured discharge against the applied 
_ pressure. The pressure-discharge relationships of the | tested emitters are Tepre- — 


sented by the specific values of aandbinTablel, 
The results presented i in Table | demonstrate the differences ii in the discharge 


= 


For the helical long- -path emitters, the value of was in the range 
- of 0.86-0. 66, depending on the nominal emitter discharge and therefore on the 
- cross-sectional area of the flow path. The value of 6b for the labyrinth-type 

emitters was close to 0. while that the vortex- was close 


00.43, 


| 

| 

| 

&g en 
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EMITTER TTER DISCHARGE SENSITIVITY 


The« experimental value of b fora given emitter determines the maximal =e ar 
in line pressure allowed in ordér to comply with the permissible 10% variation 

in emitter discharge in the field. If Q, and Q, are the maximal and minimal 
emitter discharge, then the ere: condition must be satisfied: Fo ad ye 

-— 


If we define AH = H, — H,, then Eq. 2 can be rewritten as follows: 


% 


The maximum allowable ratio, AH/H,, for the various emitters tested in the _ 


1- Values of a and b for Various Emitters" 


Helical; 1-L/h discharge 
Helical; 2-L/h discharge 


Correlation coefficient 
| 


Helical; 4L/h discharge 0.6869 


Helical; 8-L/h discharge 0. 6654 | 
Labyrinth; 4-L/h discharge a 0. 5286 
Labyrinth; 2-L/h discharge 0.5104 
Labyrinth; 4-L/h discharge | 0.4951 
Vortex; 7.5-L/h discharge | 0.4267 | 


Helical; 1-L/h 
Helical; 2-L/h 
Helical; 4-L/h discharge 
 Hidlical: 8-L/h discharge 


The. presented in Table 2 show that | for a 

10% variability in emitter discharge results in different permissible line pressure > 

‘d variability for the three types of emitters. A 15% variation in line pressure 
could be designed for the helical long-path emitters, a 20% variability for the : 


' labyrinth-type emitters, and a 25% variability in line pressure for the vortex- type 


— 
— 

_ | 
| 

R933 
13, | 0.9990 

14672 | 0.9988 

13714 | 0999 § 
2063 | 0988 

‘TABLE 2.—Maximum AH/H, for Various Emitters, as a percentage d 

AH/ 2 for Various Emitters, as percentag q 

J 
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DiscHarGe To TEMPERATURE 
_ The experimental results of emitter discharge at the various temperatures 
7 at a given line pressure suggest that the discharge temperature relationships 
for the tested emitters are essentially linear and could be described by the 


7 ( in which m and n = constants; and T= 


centigrade. The constants m and n, 


: | line pr pressure and with the cross- s-sectional a area of the flow path for the helical _ 9 
long path emitters. For the labyrinth-type emitters, a relatively small discharge | 


0.0116 0.9947 


ie. 8.8801 


2 = Labyrinth 4-L/h discharge 
48769 


—0.0222 


0.0213 


wi, 


— 


| | 
| 
| 
7 
hier lemperature, in degrees 
a computed by a linear regression analysis _ 
3 forthe various emitters. 
7 rar, 
| 
discharge 
| 


EMITTER DISCHARGE § SENSITIVITY 


sensitivity to temperature which increased with line. er was measured | 
for the Labyrinth 4-L/h emitter. discharge of the Labyrinth 2-L/h 


- Labyrinth 4-L/h emitters have shown practically no sensitivity to water tempera- 

_ ture. In contrast, the discharge of the vortex-type emitter decreased with an 
increase in temperature. Ths measured sensitivity decreased with 

The relative « emitter discharge, Qrr, asa result of changes in water temperature 
| 

“. in which Q,= = ‘the emitter discharge at a given temperature, T; and Q,= 


= the emitter discharge at a reference water temperature of a? cr” F). 


RT 


FS 


| 


= 
2 
2 
= 
= 
3 
ad 
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—Influence of Temperature on Relative ve Discharge for or Representative Emitters 


dependence of Q,, on water for representatives of the three 
types of emitters at a line pressure of 100 kPa is demonstrated in Fig. 1. Also 
_ includedi in this figure is the predicted effect of f viscosity changes with temperature 
‘The changes in the relative discharge due to temperature, Q ar? for the Helical - 
emitter were significantly less than those predicted by the ‘dependence 
of the viscosity on temperature. The sensitivity of Q ,, to temperature is expected _ | 
a me .. equal the dependence of viscosity on temperature, as assumed by Keller J 
and Karmeli (2), only when the flow is purely viscous. The writers (4) concluded % 
g the flow regime in helical long- path emitters is influenced by inertial as 
as viscous forces. Therefore, ‘the reduced of the « of 


helical long-path emitters to temperature measured i in the | 

the conclusions drawn by the writers (4), 


— The measured dependence of Q ,, on water temperature for the labyrinth-type 
emitters was small or “insignificant. These results were expected, since > the flow - 


- for the coefficient b in these emitters. The results showing a negative sensitivity — 
of Q,, to temperature for the vortex-type emitters were unexpected. The 
- measured value of the coefficient b for these emitters was 0. 0. 43. - No —— 
can presently be be given for these low val values. 


Comaineo Errects oF Pressure Temperature on Emrrter DiscHarce AN 5 * 
_ Results published by Gilad et al. (1) are used. here to demonstrate the combined = 
effects of pressure and temperature on emitter discharge in the field. Gilad» 
et al. measured the discharge and water ter temperature for each emitter in a level 
-70-m (229-ft) long plastic irrigation line exposed to the sun. Helical 2-L/h emitters 
‘spaced 1.0 m (3.3 ft) apart were used in these tests. The irrigation line had 
an inside diameter of 0.93 cm, and the line pressure at the entrance to the 
_ The line pressure at the position of each emitter, anteen’ — head losses, 
was calculated in steps using the Hazen-Williams equation, in which the exponents 
were empirically fitted for small diameter drip lines. The discharge of each 7) 
emitter along the line as a result of line pressure was computed using Eq. 
1 with the constants a and b characteristic of the Helical 2-L/h emitter. The 
_ computed results were expressed on a relative basis. The — emitter sense a 
due to to was defined as follows: 
in which Q ,,, = the discharge of emitter i tter i along t the line; and Qn, = = the discharge 
ofthe firstemitterontheline, = | 
‘The water temperature at the exit of each emitter measured by Gilad et | 
mt) was used to compute the emitter discharge due to temperature, Or. Eq. 


the constants m and n, characteristic of the Netafim 2 emitter 
‘used for that purpose. The computed results were expressed on a relative basis, 
the relative emitter discharge due to temperature follows: 


in which pets = the discharge due to temperature e of emitter i 


the discharge due to temperature of the first emitter along the line. th 
_ ‘The relative emitter discharge due to the combined effects of pressure e and 


temperature, Q,,,,, was computed using the following equation: = 

~ “RH 


The computed values of Q 2,,,Q gr, and as a function of emitter Position 
— along the irrigation line, are presented in Fig. 2. Also included in this figure 


‘ 

$ | 

| 

ey 

: 
>. 
a 
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_ EMITTER DI DISCHARGE SENSITIVITY 


is s the measured along the line published by Gilad et al. 
and expressed on a relative basis, Q py. 
The computed relative emitter discharge due to pressure, RH» ,, and the e relative 


predicted pressure dissipation along a an irrigation line. The value of Q aa 
4 gradually with distance along the first two thirds of the pipe followed a - 


ha 


itter discharge, 


vetlel 


‘lative em 


FIG. 2.—Distribution of Relative Discharge Due to Temperature, Q,,, and eee ro 
Times Temperature, Co Compared to Measured Along Irrigation 


of a given volume of water would tend to increase as it moves along a pipe 

_ heated as a result of solar radiation. Its temperature at any given time would 
be a direct function of the duration of its presence in the pipe. The variation — 
of the water temperature with distance along the pipe is expected to be inversely _ 

a related to the water flow velocity in the pipe. Since the flow velocity drops j 
sharply towards the end of the pipe, the water temperature would tend to increase 

‘The combined effects of the pressure and temperature variation along the 


resulted in a curvilinear distribution of the measured relative 
discharge, Q The values of Q em were close to the computed values of 


; _ the contrasting effects of pressure and temperature distribution on the variation | 1 


distance was more than that of Diss. The decrease i in am With distance 


ceased after about 40 m of flow. Along the last 30 m of the pipe Q ,,, ae : 


gradually, while Q,,, continued to decrease at a slower rate. The influence | 


of line pressure on Q ew Was apparent along the first third of the pipe length, 


after which Q py became sesamiae influenced by the water temperature 
predicted relative emitted discharge, computed by the 
values of Q ,,, by Op, for each emitter, described quite well the emitter discharge 
_ distribution along the 70-m pipe. The maximal variation between the measured 
and computed emitter discharge was in the range of +1.4%. __ 


Variation in water temperature along” an irrigation pipe which is ‘opposite 


in magnitude to the line pressure variation, seems to be common in irrigation 
~ pipes exposed to solar radiation. These observations suggest that emitter discharge 
sensitivity to temperature 2 as well as to pressure should be considered i in n irrigation 
system design. 

_ The discharge sensitivity to] pressure, as expressed by the. constant b, measured _ 
in the present study was found to vary from 0.43-0.86 depending on the flow 
regime and cross-sectional area for flow of each emitter type. Accordingly, — 
the maximal permissible percent variation in line pressure was computed to 
be 15% for the helical long-path emitters, 20% for the labyrinth emitters, and 


= 


25% for the vortex emitters. 


_ The discharge sensitivity to temperature was found to be positive for — 


a _ helical long-path emitters, not significant for the labyrinth emitters, and negative 
: for the v vortex emitters. In most cases, emitter discharge sensitivity to ae nae 
L The ‘measured and computed discharge distribution of helical long-path emitters. 
along a 70-m long pipe exposed to solar radiation was curvilinear. Emitter discharge 
_ followed the pattern predicted by the pressure distribution and the discharge 
7 _ Sensitivity to pressure along the first 15 m of the pipe. With an increase in 
; distance, emitter discharge became less dependent on pressure distribution and 


more sensitive to the temperature distribution and the discharge : sensitivity to ; 


‘temperature. The end result was that emitter discharge along the second half 
of the pipe increased with distance, parallel to the rise in water temperature. — 

_ Temperature variation along an irrigation pipe exposed to solar radiation is 

a result of the energy balance of the pipe and the flow velocity. The energy 

balance of a pipe is expected to change diurnally, with the season and with 

the degree of shading by the crop. Therefore, it is doubtful whether the water 

4 pans variation along an irrigation pipe could become a standard parameter 

. ‘in irrigation system design. However, based on the results from the present 


are suggested. When emitter discharge sensitivity to temperature is positive - 


work and the results cited from the literature, some general recommendations © d 


"(helical long-path emitters), the requirement for a | maximal 10% variation in 
7 discharge could be limited to the first half of the pipe. Emitter discharge is 
expected to increase along the second half of the pipe, due to a Te 


increase irrespective of line pressure variation. However, when the discharge 
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EMITTER DISCHARGE GE SENSITIVITY 

4a sensitivity to ‘temperature is negative (vortex emitters), the decrease in emitter 
_ discharge with distance is expected to be accentuated as a result of the cin emit 
- rise, especially towards the end of the line. This calls for the design of a smaller 

pressure \ variation than that normally Tequired. For -type emitters, 


variation pipe could be i in the 


‘Gilad, D., and Tirceh, G., Equipment Tests,”’ No. 81268, ICWA, 
2. Keller, J., and Karmeli, D., ‘‘Trickle Irrigation Design,”’ Rain Bird Sprinkler Manufac- 
3. Parchomchuk, J., “Water Temperature Effects on Emitter Discharge Rates,” Transac- 
tions, ASAE, Vol. 19, 1996, pp. 690-692. 
“a Tal, S., and Zur, B., “The Flow Regime in Helical Long-Path Emitters,” Journal — 
“3 of the Irrigation and Drainage Division, ASCE, Vol. 106, No. IRI, Proc. Paper 15291, 


Appennix Il.—Notation 


H = pressure; 


relative due to pressure; wi 


_ relative due to pressure and temperature; ore ee 
= relative due to to temperature; and 
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RESERVED WATER RIGHTS 
By John W. Bird,’ F. ASCE 
ITAL 
_ Because the Federal government owns large tracts ot ra in the western 
states it has a great potential interest in the waters available to the development 5 
_ of those lands. Expansion of Federal activity in those states has led to increasing — 
: Federal needs and claims for v water since the turn of the century. This has 
i to the development of a concept of a “reserved” right, the taking of 
_ state- -created water rights and a derogation of state law. _ Any agency such as 
. a municipality or an irrigation district that has a need for water and a begat 
- ‘Tight to that water must carefully consider the unwritten reserved claims | of — 
the Federal government and the limitations of those claims (if they do exist). 
‘Fer states following the riparian doctrine, the possibility ex exists that there | may 
_be existing Federal reserved rights to a quantity of flow rather than the 9 
- share ina flowing stream that is the usual rule for riparian users. In appropriation 
states a water user must consider the possibility that an existing senior — 
right may become junior to a Federal reserved right when that right is admitted 
Since there is no loss of the reserved right for nonuse or because of state _ 
laws, that right may | be valid, though hidden, far into the future. As a result, 
there have been a large number of suits concerning water rights in recent years 
that have been based on Winters’ Doctrine or the Doctrine of Reserved Rights. ’ 
It is common, in these ‘suits, to have a state either as a principle or as an 
_ interested party that is contending against . the Federal government’s claims. 
For example, there is a suit, on appeal, of United States of America v. Ti ruckee 
_ Carson Irrigation, State of Nevada et al., which is anes on rights perceived. 
to originate under the Winters Doctrine, 
"Assoc. Prof., Civ. Engrg. Dept., Univ. of Nevada, Reno, 89557. 
7 _ Note.—Discussion open until August 1, 1981. To extend the closing date one month, 
- _ a written request must be filed with the Manager of Technical and Professional Publications, 
_ ASCE. Manuscript was submitted for review for possible publication on June 25, 1980. — 
This paper is part of the Journal of the Irrigation and Drainage Division, Proceedings § 


of the American Society of Civil Engineers, ©ASCE, re a = rs March, 1981. 
ISSN 0044-7978/81/0001-0011/$01.00. 
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2S of Stee Com reserved rights. Concerning ' water rights, this is g and as such ba 
- developed from Article I, Section 8, ‘‘Congress shall have the power to regulate ; 
a commerce with foreign nations, and among the several states, and with Indian ~~ a 
_ Tribes,’’ and Article IV, Section 3 (2), ‘“‘The Congress shall have power to BY 
_ dispose of and make all needful Rules and Regulations respecting the Territory 
other Property belonging to the United States, . ... These are, 
commerce clause and the property clause of the Constitution. 
_ The first Federal suits regarding water were related to control over navigable r! 
_ waters that was reserved to the Federal government. The United States Supreme . be 
Court held that in the grant of power to Congress to regulate commerce among 
the several states, as given by the United States Constitution, the federal © 
has jurisdiction over navigable streams as far as may be necessary 
for commercial purposes (8). The Supreme Court expanded the concept of 
navigability from the thought that ‘‘a river which is navigable in fact is navigable 
in law,”’ (5) to the concept that ‘ ‘in determining whether streams are Public, 


navigable waters, the test is whether they are, or are capable of — “ 


highways between public places” (4). wm 
2 The United States made the Louisiana Purchase in 1803, the Oregon Country s 


‘Treaty i in 1846, and signed the Treaty of Guadalupe Hidalgo in 1848. The acquired _ 
lands were the property of the United States and were subsequently subdivided 
ime territories that ultimately became states with all the rights | ane ‘privileges ’ 


the thirteen original states, (62) in 1862, most ofthis > 


‘though ‘there had been a 4 steady westward migration for a number of ‘years. 

; _ None of the first migrants could claim a riparian right to needed water since 
they often were not landowners and could not be until after 1862. As a result, ‘“— 4 
_and for some other reasons, the doctrine of prior appropriation of water was ma 

perso and adopted in some areas. _ The western courts usually held that, 

while recognizing and applying the doctrine of appropriation, they were » only = 

settling rights as between individuals and that Federal ownership of the lands ad 
oe absolute vesting of such rights (11). There was a general judicial at 

_ Fecognition that the Federal government, as landowner, had a superior sight 

to the waters that were appropriated by the settlers and that there | was a 

Federal claims either by the Federal government or : by private “grantees wih 
associated superior rights. The riparian patentees could claim riparian rights _ 


based on common law. During the period that territories, states, and laws 
concerning water were being formed, Congress and the 


a agricultural, manufacturing or other purposes, have vested and accrued, 
and the same are recognized and acknowledged by the local customs, — 
& and the decisions of the courts, ‘the Possessors and owners of such . 


‘aa 


RESERVED WATER RIGHTS 


and, in in 1870, Congress clarified | (33) the situation when it said st ait anal 
all patents granted on pre- emption or homesteads allowed shall be 
‘tt ‘subject to any vested and accrued water rights, or rights to ditches and — 
and reservoirs used in conjunction with such water-right, as may have been | 


acquired under or recognized by the preceding section = ~~ 
The new states interpreted the Acts of Congress in different ways, with the 
= that they: (1) Continued with the riparian doctrine; (2) developed the 


- appropriation n doctrine; or (3) used the ‘‘California Doctrine” with 1 riparian rights" 


= 


on private lands and appropriations from the public domain forming a dual 
. _ California promulgated the theory (10) that prior to the creation of the bagi ¥ 
- states the Federal government had proprietory rights to and sovereignty over 
_ the western lands. While sovereignty passed to the states as ‘they: were formed = 
- the Federal government continued its proprietorship of the public domain. The = 
_ water rights of the Federal government were to be determined by state law 2 
: and through adoption of the common law, including the riparian system of 
_ water rights; riparian rights were created in the Federal government as a ; 
- landowner. The he Oregon courts (7) agreed with the California courts that the 
_ Federal government had proprietorial rights in the waters on the public domain - 
and that these rights survived the creation of the western states. However, 
7 Oregon viewed the Desert Land Act as severing the water rights from the : 
. public lands so that subsequent Federal land grants did not convey water rights % 
- The act was considered an exercise of the dispositive power of the Federal — 
government to the extent that it dedicated the waters to public use by a 
tion. Colorado viewed the sovereignty of the Federal government as separate 
. § from proprietary rights. Thus the transfer of sovereignty to the states as they 
_ were formed transferred all the power to control the waters. The courts interpreted — 
— the Act of 1866 along with subsequent pag enactments to give the —— 
consenttothisinterpretation, j= 
Since Congress did not adopt any Federal the for the disposition or 1 determina- _ 
4 tion of Federal water rights, the newly formed territories and states took the 
_ function of establishing water rights based on common law, statutes, and court 
- decisions, and in doing so they were exercising their traditional jurisdiction 
over water. Each appropriation state has similar basic legal concepts: (1) Beneficial 
use is the measure of the existence and scope of the right; (2) ownership of 
the land is not considered a basis” for a water right; (3) the right may, but 
need not necessarily, be appurtenant to the land; (4) the appropriated water 
_ may be applied at any place where it is needed, regardless of the distance _ 
4 from the source; (5) diversions are protected; (6) the sights of the prior appropriator ‘4 


the ov owner of | a senior can take no more water than was "was necessary 
_ for his original need; and (9) the right extends only to the use of the water, 
there is no right to the corpus of the water. Apparently no one at the state 
level considered the possibility of a Federal sovereign right or how such a — 
= might bi be established. 1 The Act of of 1866 wi was considered to be a ahaa 


| 
| 
domain, 


a In pest, , the D Desert Land Ac Act of 1877 7 (34) sa said 


_ . the water of “all lakes, rivers and | other sources of water idle 
= the public land and not navigable, shall remain and be held free + 
- for the appropriation and use of the public for ‘irrigation, 1 mining, and 
manufacturing purposes, subject to existing rights. . . 
which had been interpreted to be a dedication to the “public of all 
i ia the waters of the public domain (7). This interpretation separated water from 
7 mane & and was apparently upheld later when the Supreme Court held (2) that, 
7 “2 |, all non- navigable waters then a part of the public domain became publici 
ji, [of public right], , subject to the plenary control of the > designated states 
rf ’* Using these interpretations the various states passed laws and developed | 


- 2 will be recalled that the first suits concerning 1g water that included the Federal 
government were in relation to defining and determining navigable waters. This 
was changed and expanded in 1899 when the United States sued to restrain 
the building of a dam across the upper Rio Grande River based on the grounds ; 
that navigation would be impaired by that construction because of resulting 
_ diversions from the Rio Grande (23). It was the holding of the Supreme Court 
that, although the states have the right to adopt riparian or appropriation ned 
and regulations as to flowing waters for such purposes as they deem proper, 
4 limitation must be recognized in that in the absence of specific authority _ 
from Congress, a state cannot by its legislation destroy the right of the United 
States, as owner of lands bordering on a stream, to the continued flow of 
its waters. It is well recognized that the superior power of the Federal government 
to secure the uninterrupted navigability of all navigable streams within the limits a 
of the United States. Further, that power may be used to prohibit any obstruction _ 
to the navigable waters of the United States. ‘This power v was s and is reserved 


= the non- navigable portion of the Rio Grande, the Federal government had re 
_ the authority and power to extend its control from the navigable portion to 
_ the non-navigable headwaters in order to protect its sovereign rights downstream. — 
- While it was contended that the Act of 1866 and the Act of 1877 represented 
a Congressional grant of complete authority © over non- navigable waters to the 


4 

: 


- gelease that power (23). Further, even if the acts had the contended effect, 
the subsequent Act of 1890 (30) prohibiting obstruction to navigation must be 


7 that the navigation clause will strengthen federal authority for the project. aa 
_ This decision that there were reserved powers in the Federal government 
over waters that were in fact non- aavignate was extended and applied in _ 
Winters case. = 


: 
: 
| 
{ 
4 
1€1G tO since there is no limitation On sudsequent Congressional actions 
_ regarding navigation. One result of this ruling is that Congress will often claim : 


RESERVED WATER RIGHTS RIGHTS 
. ae series of cases involving the water rights of Federally owned /administered 
i Indian reservations illustrates Federal claims based on Federal proprietary q 
_ interests and implied reservations. Indian reservations were created by treaty, 
proclamation, or executive order while the west was still largely federal territory. 
_ Reservations were typically established to change Indians from a nomadic to 
an agrarian people (26). After formation of western territories and states the 
available waters were usually appropriated by non-Indian settlers, leaving little 
water for use on the reservation. Other reservations, created for other purposes, 


also had little or no water available for their use after the available surface — 
water had been appropriated for other existing needs. __ 
To protect its program of Indian agrarian development the Federal government | 
asserted a claim of “‘reserved’’ water rights against non-Indian appropriators 
beginning in 1908 with Winters’ in Montana (27). In the Winters case, the United > 
States claimed for the Indians living on the Fort Belknap Indian Reservation 
: on the Milk River that, when the United States created the reservation, by 
- - implication, it also reserved sufficient water to carry out the purposes of the 
reservation a right to water for the purposes of irrigating their land, even though "i 
the treaty that reserved land for them made no mention of water. In 1899 . 
water was diverted from the Milk River by Winters and other defendants via 
dams and canals. As a result of the conflict between settlers and Indians the 
United States Supreme Court held there was an implied reservation of the amount 
- of water necessary for successful cultivation of the Indian lands. The reserved 
r = had a date of priority of the formation of the reservation and was superior — 


to later appropriations made under state law. The ‘Supreme Court held that 


& the Desert Land Acts of 1866 and 1877 were “4 ee to the lands or 
waters of Indian reservations, 
‘The reserved-rights doctrine applies to all toftes: reservations, even though > 
Di many of them were created after 1866 or 1877 and regardless of how the reservation 
in question was: created. The source of the reserved waters includes waters - 

Sane upon, flowing through, or bordering Indian reservations, with the date 7 

¢ reservation of water coinciding with the date of creation of the reservation. 

The question of how the Federal water rights could be released from public 
4 lands that were later turned into reservations with implied reservations of water 7 

“ The Federal courts have based their decisions concerning Indian water ' rights 

on one of several theories. The theories used are: (1) The United States reserved 

- (implied) water rights in the treaty that created the reservation (1); and (2) 

_ the Federal government reserved the land and water (implied) from the public — P } 
domain for the Indians (25) (no treaty). It is unclear whether the courts ine 

- the latter cases regarded the reservations created by | the Federal vet likes 

as made in the exercise of the property power or the treaty power, but either 
Reserved wat 
withdrawn from the public domain. The Public domain those lands 
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owned by the United § States that are subject to entry, dunenth: or private 
appropriation under public land laws, while reserved or withdrawn lands are 
taken from the anude lands a and are je not subject to to entry, dieposal or peivate j 
and i is not a public right but is s considered to be a . private e right held in trust 
by the Federal government for the benefit of the Indians. When the reserved 
right was first developed it was considered by many peagie to be an incident 
_ associated with one suit and then, after several suits, with Indian reservations 
only. This has changed as a result of suits concerning reserved lands that have 
been determined to have reserved water rights associated with them (12,17). 
- These other reserved water rights are public in nature, but all of the reserved 
_ have a priority date based upon the date of creation of the reservation. 
The Indians’ reserved water right, when used for irrigation, appears to ie 


appurtenant to the land, which is similar to state-created water rights. Those 
reserved water rights cannot be > lost | by nonuse under state laws, nor by legal 


Bs through condemnation or state statutory enactment, or private appropria- — 
- tion. While the Indians may be directed to follow state laws in some cases, 
that statutory language should not be construed as an abandonment of prior 
existing rights by the Indians unless that intent was clearly expressed in the 
statute. The statutes only provide for procedural filing under state law, but — 
do not limit the reserved right. The overriding power of the federal sovereign — 
7 under the supremacy clause (Art. VI, Cl. 2) of the Constitution is the source 
; of the protection of Indian water rights from encroachment. Unless Congress 


3 clearly consents to state jurisdiction, the right cannot be set aside, overridden, . 

denied, or otherwise be affected by state substantive law (1). The decision 

that the implied reservation may be for future unidentified use (1) constituted 
a significant departure from appropriative water law that has resulted in —- 
erable opposition from the states and non-Indian water users as there is no 
well-defined measure of the amount of reserved water for Indian uses = 
so non-Indians have no assurance of the quality of water left for their use. __ 
From the first cases the water reserved for use by the Indians was for agricultural 
purposes, uses, and needs. In Arizona v. . California this was modified to be 
measured in terms of the ‘‘practicably irrigable acreage’’ on the five reservations 7 

_ involved and not just for the land irrigated at that time. Reservations of water 

3 were determined to have been created for the purpose of enabling Indians to 

‘ develop a viable agricultural economy and other uses, such as those for industry, 
which might consume substantially more water than agriculture, were = 
contemplated at the time of the creation of the reservations. Indeed, in the 

_ middle 1800s a common aspiration in the United States was to have one’s own * 

_ land to farm. It was made clear that when an Indian reservation is established 

_ to provide an Indian agricultural economy, the measure of the associated reserved 
water right will include that amount of water necessary to irrigate the practicably 
irrigable acreage and to satisfy related uses. . The question as to whether the 7: 


2 answered. Furthermore, nothing has been said to date about an Indian reservation 
which may have not been created to establish an agricultural economy. Does © 
a creation of a reservation that may have, as one of its original purposes, 

the preservation to the Indians’ benefit | of the fish and waters of a lake, have — 


| 
| 
| 
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reservation of to maintain tet lake elevation | 
fishery? This question is currently being litigated (24). 
3 Two standards have been promulgated for determining the uses that have — 
4 1. Those uses necessary to fulfill the purposes contemplated when the 
4 Teservation was created (used in Arizona v. California and Winters’ Rae : 
a All possible uses, s, including those that may appear in the future without — 
reference to the purposes for which the reservation was created [United cect 
v. Ahtanum Irrigation District(14)}. 
The first method permits a specific amount of water to be identified and yt 
from other users since there would be an immediate quantification of the reserved | 
right. The Indian’s rights remain uncertain with the second method and, to 
- some extent, unprotected since the courts could conceivably require the Indians 
- totake compensation for their claim if unappropriated water is no longer available. 
_ The suits concerning waters for reservations, for power, and for other purposes’ 
lead to the McCarran Amendment to the Justice Appropriation Act of 1952 
_ where the sovereign immunity of the United States was waived (35): _ oh well ah 
any suit where it ‘appears that the United States is the 


This. has S permitted some states to bring suit to include the { United States, through 
its a agencies, in state court - water-rights ; adjudications to disclose and quantify 
any reserved rights where jurisdiction can be established under the etna 
_ In retrospect it may appear obvious that if Indian reservations had reserved 
- water rights, then all Federal reservations should have some implied Federal 
reserved water right. In practice, it appears t that almost ev everyone considered 
the: reserved right to be confined to Indian 1 reservations and to be, in fact, — 
a strange quirk of law associated with Indian treaties or Indian reservations, 


or both, and not to be associated with any non-Indian reservations. 
y 


difference withdrawn, ar and public lands was emphasized 
_ when the Supreme Court upheld the authority of the Federal Power Commission - 
a to license construction of a private hydroelectric project on Federal reserved — 
lands in Oregon despite state objections over the ; Pelton Dam (6). Since the © 
Stream involved was non- -navigable, there was no “Ro assertion of the navigable 
servitude that is reserved to the Federal g« government. The Pelton Dam case — 
; considered to be a landmark by some because the Court concluded that 
_ superior Federal authority, under the power of Congress to regulate the use > 
and disposition of Federal property, existed to license the power project in 
violation of state law. The view of a lower court that the Desert Land Act > 
constituted an irrevocable control over non-navigable waters to the state was 


_ implicitly rejected with the statement that the Desert Land Act does not apply 


| 

- 

1 


power purposes—well af ter ‘the Desert Land Act enactment. It was not explained — 
fc how the creation of a reservation withdrew the property reserved from the _ 

operation of the Act. In the dissent it was stated that while Congress could | 
A have revoked the Desert Land Act as far as reserved lands were concerned, — 


= for geaiente on non- -navigable v waters “were to be acquired by following 
state law. But this view did not se | Since the seiay Court did not fully — 


apart from ownership, to control the waters passing through its lands, subject 
to payment of compensation for interference with existing rights. Since the 
Desert Land Act severed the waters from the lands so that ownership of the 
_ waters did not pass to patentees of Federal lands but remained Federal property 
subject to state control under state law, the combination of Federal ownership A 
of both reserved lands and unappropriated waters may have controlled the Court’s __ 
decision. If this power of control exists, it must extend beyond the boundaries _ 
Ss the reservation to insure that waters required for the reservation continue a 
to flow down to it [see Rio Grande (23)]. It should be noted that this view | 
would appear to be contrary to the holding in Kansas v. Colorado (9) that 
the property power gives control only over Federal property and not unappro-- 


_ ‘The uncertainty concerning the reasoning behind the Pelton decision raised 7 


the west. This was augmented by the Hawthorne Naval Reservation case (12) 
_ where the Federal court held that in maintaining a national defense reservation, — 
a valid Federal purpose supported by both war and property powers, the United 


control ov over r all non- n-navigable waters flowing through Federal reservations in 


~ complying with state law. This was an application of the j property power reserved — 
* that was extended to subsurface water, + Nga 
In 1956, in the State of Washington (14), it appeared that the Indian’s water 
“rights were extended when the court said that Indian rights were not limited 
“tothe the use at any given date, but rat rather 
extended to the ultimate of the Indians as those needs require- 


_ The court then upheld an agreement inenaen. wii white settlers and the erat 
of the Interior, executed in 1908, that limited the settlers to 75% of the Ahtanum 
waters. The water right was again expressed for irrigation purposes and implicitly — 


_ Since the construction of the Colorado Aqueduct, California has been diligent 
4 in appropriating water from the Colorado River. In an effort to determine and 
protect its water rights i in the Colorado River, Arizona filed a suit against Calif California 


Court mav have recarded nower nronerty as the sovernment nioht J 
[. 
: 4 
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in 1952. Nevada, Mexico, Utah, and the United States were added as 
= = to determine how much water each state has a legal right to use from 
the Colorado River and its tributaries. Simon H. Rifkind was appointed = 
_ Master and reported his findings to the Supreme Court in 1961. In those 
"proceedings the United States asserted claims to waters in the main river and 
_ in some of the tributaries not only for use on Indian reservations, but also 
for us use in national forests, recreational and wildlife oenren other government 


and law that 


of the reserved lands (1) or to supply other primary needs. Arizona argued 
that: (1) The United States had no power to make a reservation of navigable — 
waters after Arizona became a state (1912); (2) that navigable waters could 
ae be reserved by Executive Orders (some of the reservations were created _ 
: by Executive Order); (3) that the United States did not intend to reserve water 
7 for Indian reservations; (4) that the amount of water reserved should be measured > 
by the reasonably foreseeable needs of the Indian living on the = 
4 rather than by the number of irrigable acres; and (5) that the judicial doctrine 
be used to divide the water between the 


Supreme Court answered Arizona stating that equitable 
7 is used to solve disputes between states and an Indian reservation that is not 
7 a state. The argument that the United States did not have the power to reserve 
navigable waters was based on cases that concerned only the title to lands, 
on ‘not waters, adjoining or underlying navigable waters. The power of the United 
States to regulate navigable waters is not reduced or limited. The arguments 
regarding the intention to reserve water and the amount reserved were rejected 
with reference to Winters (26) with the statement that the water rights were _ ; 
reserved effective as of the time the reservations were created. 
-. For the first time the Supreme | Court extended the reserved doctrine of water 
“rights from Indian reservations to all reserved or withdrawn lands such as wildlife 
— and national forests. Since the Court considered only claims on behalf 
_ Of Indian reservations, it is necessary to refer to the report of the Special 
Master to determine the basis for extending the reserved rights doctrine to 
_ other reservations. The Special Master first determined that the United ‘States 7 
had the power to reserve water to fulfill its purposes in cr 
reservations involved in the suit; then he determined that the | circumstances 
‘surrounding their reservation demonstrated the intent of the United States to 
reserve water for use on the reservations involved. These reserved water rights 
_ are subject to water rights that have been previously acquired pursuant to state 
law, i.e., that have a priority older than the « date the reservation was created. — 
President H. S. Truman reserved 40 acres of land as a national monument, 
in 1952, under the Act for the Preservation of American Antiquities (13) which | 


authorizes the President to proclaim as national monuments “‘objects of historic 


_ or scientific interest.” That 40 acres included Devil’s Hole, a deep limestone 


cavern including an underground pool inhabited by a unique species of desert 


— 4 when the United States created those reservations, or added to them, it reserve . a 
only land hunt aleo the nce of enonoh water to irrigate the irrigahle 
a 
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for the sails is a rock ledge in the pool and if the grouad- -water level drops 
_ Significantly that rock shelf becomes exposed to air with a resultant decrease 
in spawning area and the p potential extinction of the pupfish. As a part of the 
- preamble to the proclamation proclaiming Devil’s Hole as a national monument, — 
President Truman said, ‘“Whereas the said pool is of cummnating scientific 
_ importance that it should be given special importance. . 
In 1968, Cappaert, a local rancher, began pumping ground-water from wells 
located | about 2-1/2 miles (4 km) from Devil’s Hole. The pumping pn 
during the irrigation season between March and October of each year and Cappaert 
was the first to appropriate ground water in the area. The National Park Service. 
filed a protest to his ground-water application requesting that it be denied until 
- it was determined which wells, when pumped, were causing a serious decline 
-@ water elevation in the pool. The state engineer declined to postpone i 
decision because ‘there was no recorded Federal water right with respect to 
Devil’s Hole. In a suit in Federal court the National Park Service claimed + 
that the establishment of Devil’s Hole as a National Monument reserved the 
unappropriated waters appurtenant to the land to the extent necessary for the 
_ purposes and requirements of the reservation. Since neither Cappaert nor anyone 
else had perfected water rights as of the date of the reservation, the _ reserved 
“ rights had first priority. An injunction was issued in Federal District Court 7 
; i in 1973 limiting the pumping from designated wells to control the pool elevation 
to an acceptable level (17). On appeal, the Ninth Circuit Court of Appeals 
affirmed the District Court decision (17) and held ‘‘the United States is not 
bound by state water laws when it reserves land ie the public domain.’ - 
The United States Supreme Court affirmed that decision in 1976 when it considered _ 
the water in the pool to be surface waters and not ground water; as a result, 
_ there may be a little doubt on whether or ‘not the reserved right does extend. 


i In 1971 the Supreme ‘Court declared (20) that ‘ “it is clear from our cases 
that the United States often has reserved water rights based on a 
from the public domain.’ As we saw in Arizona v. California, 373 U.S. 546, 
the Federal government had the authority both before and after a state is admitted 
into the Union ‘‘to reserve waters for the use and benefit of federally reserved — 
lands.’’ The federally reserved lands include any federal enclave. From the : 
Eagle County (20) affirmation of Arizona v. California it is clear that water 
rights have been impliedly reserved to serve not only Indian reservations but 
also any Federal enclave cr created by -Teserving or withdrawing lands from the 
_ Perhaps the principal question that arises is—what reserved lands? The 
following i is a partial list of some of the reservations; it does not include Indian | 
‘reservations, power site withdrawals or military installations since these have 
“been previously mentioned. Congress specifically authorized the reservations a 
s some public springs and water holes in 1916 in the Pickett Act (36) where — a 
‘lands containing waterholes or other bodies of water needed or used by | the 
‘public for watering purposes... reserved . for public purposes. 


‘ 
| 

| 


But Buress of Land Management is prohibited making wit withdrawals 
under the Pickett Act by its Organic Act (40), 
_ Many springs and water holes were reserved by an executive order for Public 
. 107 in 1926. While numerous other specific 
_ were made both prior and subsequent to y the 1926 order, these reserves were 
generally local in character or otherwise minor. This withdrawal order =a 7 
_ designed to preserve for general public use and benefit unreserved public lands’ 

a containing water holes or other bodies of water needed or used by the public 
for watering purposes. This withdrawal reserved both land and water for public _ 
use. The quantity reserved was the total yield of each source to be used where 

: } it was needed for domestic, livestock | watering, irrigation, soil, and fire -and 

_ erosion control. While the 1926 order does not interfere with vested rights under 
State law prior to the 1926 withdrawal date, if a spring or water hole is ell 
after 1926 it does apply. Similarily, if a state water right was not acquired — 

Bes after 1926, that water right i is ineffective against the 1926 withdrawal. —- 


. . the boii of oil shale, and lands « containing such deposits pa 
be the United States, be, and the same are hereby, temporarily cs 


purpose” that was Mexico ( (21) it 
appears that only enough water was reserved as was reasonably necessary for 
the “‘purposes of investigation, examination and — and not for actual © 
‘Taylor Grazing Act (37) directed the of Land Management to 
_ manage the public domain for grazing purposes. It did not ‘Teserve any land | 
~ from the public domain so, presumably, no reserved water rights were created 
by the act. Since the Classification and Multiple Use Act of 1964 (39) was 
to be consistent with and supplemental to the Tayeee Act, it should create 
4 ~ National parks created before 1916 were established by legislation that at usually 


shall mé make r providing for. the preservation from injury 


In 1916 the National Park Service Organic Act was with | Pry 
_ purpose “*. . . to conserve the scenery and the natural and historic objects — 
and the wildlife therein . . .”” Whether the national park was created before 
or after 1916, water is reserved for the purpose of the park—protecting forest 
"stream flows, “lake” irrigation, wildlife “and bind watering, po 
campground uses and maintenance, public facility uses, etc. (21). If additional — 
water is needed, or if water is needed for secondary purposes, the ied 
‘Pet Service may acquire additional water a by f following state law. yeay 
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q 
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— shall not be construed as a reservation of the waters of such 

for purposes other than those specified . . . or in quantities greater than 
mecessary to accomplish these purposes. = 


It may be inferred that the amount of unappropriated waters that are reserved 


relating to the control, use, or distribution water used 


a This has been ain ii interpreted by the United States Supreme Court (3) 


iy general, the implied Federal ‘reserved right is ‘for enough w: water to ete 
the specific purpose for which the land is reserved. Nonreserved rights are 

_ those for secondary uses, or for congressionally authorized management objec- 
tives on multiple-use public lands. In the case of nonreserved | rights the a 

government would follow state law in acquiring a water right. 

g As a rule, reserved water ‘rights were adjudicated in Federal courts until. 
the passage of the McCarran Amendment in 1952 (35) which granted a limited — 

} waiver of Federal sovereign immunity to permit suits against the United States 
in stream adjudication where the United States’ water rights are involved. The 
Supreme Court interpreted (20) the amendment to include the Federal water - 
rights of federal I non-Indian r reservations and enclaves, as well as Indian r reserved — 

j rights (15). From these cases it is now possible for the states to initiate adjudication 7 
for determination, extent, existence, scope, and measure of reserved rights in : 


the state courts without having to wait for the Federal government to initiate 
“Such a suit, 


- by Indian reservations for nonagricultural uses. There are, however, these aiiee A 
a cases that demonstrate that the measure may be that amount of water necessary — 
to fulfill the particular purpose for which the water was impliedly reserved. 
First, where a water right is asserted for the purpose of sustaining a viable 
.- in a desert lake and its supporting stream, the United States claims a 
sufficient water to maintain the present level of the lake over a long period 
of time, along with sufficient stream flows to sustain spawning runs and to — 
preserve the instream habitat for the fish and their fingerlings (24). This case 
‘involves Pyramid Lake, a large desert lake enclosed entirely within the Pyramid 
te Indian Reservation. Second, where coal mines exist on an Indian reservation, a 
3 claim is for sufficient water to bring the c © a marketable state on 


— 

22, 

: scientific, biotic, or historic features that led to the river’s inclusion as a component 
The Reclamation Act of 1902(38) provided that 
nothing in this Act shall be construed as affecting or intending to 

4 
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the Ute, Mountain Ute, an and Southern Ute Indian reservation in southern eine 

(15). Third, if preservation of the ecology of a stream is the purpose of a 

reservation, the claim if for a minimum flow of water sufficient to maintain = 
_ the environment of the stream and its rT values | on Chamokane Creek 

on the Spokane Indian Reservation (16). _ 


for the purpose “of carrying out the congressionally authorized 
of the Department of the Interior. Examination of the Solicitor’s Opinion leads 

5 one to conclude that Federal nonreserved right is limited to differences between - 
_ Federal and state requirements for an appropriation. ‘Under state law the Federal 

government the > right to file for an appropriation of water the 


- application or actual beneficial use. These appropriations are not to adversely _ 
affect existing rights established under state law. The procedure, then, would _ 
_ follow state law. The purpose or beneficial use of the right might not follow _ 
“ state law, since the definition of beneficial use varies from state to state and : 
since some states may not | accept instream uses as beneficial. While following — 
4 state procedures to. acquire a water right, , the Solicitor’s Opinion states 
that the Federal government will not be constrained to state limitations that 
s tend to constrain Federal definitions of beneficial use on public lands. 
_ The area of conflict is in determining if this right exists at all and, if it does, — 
it will probably be limited | to different interpretations of beneficial uses, with 
appropriation rights. that. may ‘not be able to be certificated to any other =i 


Fururne se 

_ The doctrine of reserved water rights has gr grown from a concept of reserved 
‘Federal powers in both navigation rights and property rights along with its — 

“treaty- making powers. The admission of a state into the union on an ‘ “equal 
footing” did not divest the United States of its plenary control over the waters — 

_ (1), but the Supreme Court has recognized (3) that there is a conflict because 
the state, upon admission to the union, , should have ‘‘exclusive oo 

over the unappropriated waters of their streams. While it may be ‘difficult or 

: impossible to have agreement between the states rights and reserved- -rights” 


- doctrine after statehood, the reserved- rights commas! is nevertheless law at this 
_ Reserved rights always begin as implied rights that must be quantified. There 
-y still many unanswered questions concerning reserved rights such as: (1) 
Can the United States reserve water rights for acquired, , rather than withdrawn 
lands; (2) are there reserved rights for minimum stream flows; (3) does a change 
_ in the purpose of a reservation destroy or change the date of a reserved right; wa 


(4) what limits are placed on the purpose of a reservation; and (5) does a 


| 
i} 
The Solicitor's Opinion, released June 23, 19/79, refers in part to Federal 
a 

| 

«water right extend olf the reserved land to a distant stream or water 
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D DRAINAG E DIVISION” 
PARAMETER Estimation 


UsING KALMAN 


‘By Subhash Chander, N. Kapoor,” 


INTRODUCTION 
‘The of an aquifer as a potential ground-water resource 
- depends largely on the correct determination of its two inherent characteristics: 
its ability to store and its ability to transmit water. In case of leaky aquifers, 
the leakage factor of the overlying semipervious layer must also be known. 
_ These parameters are estimated by analyzing the observed head variations in 


an aquifer” test. The measured heads i in the field do not “necessarily satisfy 


system and noise in the measurements. Conventional approaches for parameter 
estimation from pumping test analysis (5,9) do not account for such uncertainties, 
_ and do not give any index of reliability of the parameters. = 


_ A computationally-efficient nonlinear filter, known as an iterated extended 


Kalman filter (IEKF) has been used to get the “parameter estimates both for ¢ 
nonleaky and leaky aquifers, along with the confidence limits, in the a) 
of modeling and observational errors. 

Nonleaky Aquifer. —Under ‘the Dupuit- -Forchheimer assumptions, the aaa 

_ differential equation representing the radial flow of groundwater i in a confined : 


' Prof. in Civ. Engrg., ‘Indian Inst. of Tech., Delhi Khas, Delhi, 
Asst. Prof. in Civ. Engrg., Indian Inst. of Delhi, Hauz Khas, Delhi, 
, Indian Inst. of Tech., Delhi, Hauz Khas, New ‘Delhi, 
a Note.—Discussion open until August 1, 1981. To extend the closing date one month, a > 
a written request must be filed with the Manager of Technical and Professional Publications, _ 
ASCE. Manuscript was submitted for review for possible publication on July 2, 1980. 
_— This paper is part of the Journal of the Irrigation and Drainage Division, Proceedings - 
the American Society of Civil ‘@ASCE, Vol. 107, IR1, March, 1981. 
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in which T = the » tranemicsivity, in square meters per day; s = the unten, 
in meters; Ae = the distance, in meters, from the pumped well; S = the storage 
‘coefficient; t = the time, in days, elapsed since pumping started; and g = 
- the rate of withdrawl per unit area, in meters per day. Eq. | assumes that 


the pumping well f layer during the - Theis’ solution 


Ww 


= avai well function for a nonleaky aquifer. Eq, 2b may not be ally representative — 
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of the grou ground- water system em under analysis, and in ‘there may 

in the observed drawdowns. Due to these likely errors the estimation of f parameters 

T and S can be carried out using the IEKF (3), as in the following, = = : 
Defining the state vector as X(t, ) = [T S] ae in which [ ]’ = the transpose 


‘The monsurement equation 


i 


in Which “measurement function h() is nonlinear in T and S; the 
total error in the modeling anc and observation; it is white gaussian with zero mean — 


and covariance defined as 


| 
FIG. 2.—Estimation Error Covariance of Parameters (Leaky Aquifer) 
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& hates 
7 which R = mea 


wt 
Starting with the initial estimates for state as 


It.) = 2) = 


and for the estimation error covariance as 


P(t, ) = P(t, )=E{[X(t,) - 
the predictor-cc  algorithen for the state estimates is given 


= = silty + K(t,.,)( (teas) A 


- in which K( ) = the filter gain or Kalman gain matrix 


| 
= the two state variables, 


performing a a local ite iteration, retaining a at the. same time, the recursive 


filter structure (1). The local iterator for Eq. 11 is given as: ys aa 


in which = sumer of i iterations. 


The iterator starts with y, = silty ). It gain K, 
function A, and the transition matrix H, on each iteration, and terminates when 
there is no significant difference between the consecutive iterates. The result — 
of the last iteration, i.e. Tei is taken for the estimate and the 
_ estimation error covariance in Eq. 12 is computed on the basis of this estimate. _ 
- Measurement noise covariance, R, is adaptively estimated iia to Sage 

and Husa (6) as: 


in which the measurement residual 


and § = the Kronecker 
| 
| 
| 
» 
Variables T and S in the present case. The 
| | 
| 
| 
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the Dupuit- assumptions, ‘the partial differential equation n for ground- 
water flow in a confined leaky aquifer of constant depth is: fe ==y = 


the teenage factor in meters. Eq. 17 assumes that the pumping 
* fully penetrates the aquifer r layer during | the test. Hantush and Jacob @) ; 


Aquifer Parameter Estimates for Nonleaky Ac Aquifer Using Different 


Square for the 
analysis per oy icien servations 
0.123 158 x 107° 
0.311713 x 107° 
0.184373 x 107° 
0.110 926 x 
ef ab 


= = = 
Aq 


Leakage 


analysis” | per day coefficient meters 


curve 1,729, | 6 18 999 x 
on 
discrete 

| 

1,680 + so | 0.00 0002 "0.32 426 x 

1,658 668 0.38215x 10 


. oy observation well is at 90 m from the pumped well for all methods except tS 
Rushton’s discrete numerical model, for which the observation wells are at 30 m, i 
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“gave the fer Eq. 17 


h = well feaction | for a leaky aquifer. Eq. 18 is alternatively 


4Tt’ 


Eq. 19 is expanded as: 


which K,( ) = the modified Be 


"presence of mn observational noise can n be solv ed by u using 
the IEKF as explained earlier (Eqs. 3-16), with the State variables and the a 


measurement function given as follows: 
State vector 43 = 


easure 


(u\ry 


‘ 
q 


In this case Li are three state variables, 


hawla (7) using Theis’, Jacob’s, Chow’s, and Kriz’ methods were selected 


test data for aquifer already analyzed by Sharma and 
parameter estimation by the IEKF. The Dalem test data ha 


q | — 
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been earlier analyzed by Krus Kruseman and de Ridder (4), using Hantush and Walton’s — 


methods, « and by Rushton an and Chan (5), using a discrete numerical model, was — 
"chosen for the leaky aquifer parameter estimates. 


_ A lithological cross section of the test site in the leaky aquifer has Holocene | 


as the overlying semipervious layer and Kedicham (lower Pleistocene) as the 


impervious basis layer. During the pumping test, the screen was open only © 
- in the Kreftenheye (upper Pleistocene) which was 11 m-19 m below the ground — 


surface. This formation is of medium coarse sand with 0%-2% clay. = od} 


_ Initial estimates for the parameters are normally prescribed on the basis of 


=. ‘the geologic information of the test site. It is felt that faster convergence may — 
(ae be achieved by choosing the initial estimates on the basis of some information 
: about the system behavior to pumping stresses. The initial estimates for the 
ae and ‘storage coefficient were found out by u using the first two 
observations for drawdown and applying Jacob’s relation (4) to the two drawdown — i 
values. The initial estimate of the leakage factor in the case of the leaky aquifer 
was Calculated by using the last available observation and applying the approximate 
_ Steady-state drawdown relationship for the leaky aquifer as used i in the Hantush ~ 
_ ‘The filter provided the parameter estimates along with the estimation error 
aa covariance, and the sum of the residual squares for the difference between 
4 the observed heads and those calculated using filtered parameter estimates. 
Different sets of data were analyzed for estimating aquifer parameters. Successive 
a values of estimation error covariance for the parameters for nonleaky and leaky 


=~ and the residual squares, as obtained by using various methods, are given in 


Table | for the nonleaky aquifer, and in Table 2 for the leaky aquifer. The ‘ 


filtered aquifer parameter estimates along with their upper and lower limits — 
for 95% confidence level Gaming aGa Gaussian distribution) are given in Table 3. a. 
_ TABLE 3.—Confidence Limits of Aquifer Parameter Estimates Using i - wi 


20) esagby; ignuwihhs dass ber al soplatoss 


square meters per 
"Storage coe coefficient 


Transmissivity, in 
square meters per 
coefficient | 0.00 168 =| 0.00174 | 0.00180 


= aquifers are shown in Figs. | and 2, respectively. The estimated aquifer parameters — 


7 


| 1 
| | 
oi 
0.000029 


Comparing the the aquifer parameter estimates using the IEKF method with 
obtained using the known techniques, the following points are observed: sh 
7 1. The sum of residual squares computed from the ‘difference of the observed © 
and the calculated heads is the least for the parameter estimates resulting from 
& IEKF method for the data segment 6-15, as compared to those obtained 
: ss Error covariance plots for the coefficients of storage and transmissibility 


for the nonleaky aquifer (see Fig. 1) show a rapid decrease for the first few 
observations, and thereafter the plots indicate a marginal decrease with 
each observation. It shows that the parameter estimates tend to converge after 
_ a few observations and additional observations do not significantly change the 
3. Also in the case of the leaky aquifer (see Fig. 2), the error covariance — 
for transmissivity and storage coefficient decreases rapidly for the first first few 
observations; after that the decrease is not as rapid. The error covariance of 
_ the leakage factor decreases only marginally for the first few observations - 
_ which may be attributed to the fact that there is some time lag between the — 
start of pumping and the leakage effects influencing the drawdown. > A 9 
From the aforementioned analysis, it is concluded that: 
a 1. The IEKF provides ‘the estimates of tl the aquifer parameters (r and S for — 
_ the nonleaky aquifer, and 7, S, and L for the leaky aquifer) in the presence 
system uncertainties and noisy measurements. 
2. The filter also provides the confidence limits for the pa parameters. — Ar 
(3 The technique is sequential, and each additional observation reduces the 
error covariance of the parameter estimates, although only a few observations 
are sufficient to provide a reasonable estimate of the aquifer parameters. __ 
4, The subjectivity involved in type- curve ; matching for the parameter estima- 


tion is eliminated using the IEKF method. | Mmevecnt? 
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‘The following symbols are use in in this paper: 
H ‘measurement transition matrix; 
measurement function; 


Kalman gain matrix; 


My modified Bessel’ s function o of second kind and zero © order; 


rate of pumping, in cubic meters per eu: oe 
rate of withdrawal per unit area, in meters per day; wy ' 7 
R = measurement noise covariance 
= distance to observation well, in meters; 


= transmissivity, in square are meters per day foe eS 


well | function for fora leaky aquifer; ween oe 
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AND WATER- TYPE EFFECTS 


Yo. on GROUND-WATER QUALITY q 


By Donald L. Suarez" and Martinus Th. van van Genuchten* 


_ Leaching to prevent the of soluble salts in the crop rootzone 


F is essential | for s sustained irrigated agriculture. Under present ir irrigation practices, — 
— however, many irrigated lands are leached more than is necessary to prevent 
_ yield reductions from excess salinity buildup. Reduced leaching should not result 
7 in reduced yields in many irrigation projects in the United States(10), 
_ When leaching fractions are reduced, the total salt load of the drainage water 
is reduced (3,5), although the salt concentration in 1 the lower part of the rootzone 2 


‘salt content of drainage water under certain circumstances. In a , previous ; study, 
"the first author and Rhoades (7) analyzed the effects of improved irrigation 
_ efficiency on ground-water salinity for closed river basins. Depending upon 
water chemistry, changes in irrigation efficiency at steady-state water and salt — 
fluxes may or may not affect downstream water quality. They classified irrigation ; 


waters into three groups: Type 1—waters initially undersaturated with CaCO; 
_ Type 2—waters initially saturated with CaCO,; and Type 3— waters nearing» 
- ‘saturation with gypsum and saturated with CaCO,. Improved irrigation manage- _ 
ment resulted only in a slight reduction in downstream salinity if Type | water — 
was used, in no steady-state reductions in downstream salinity with Type 2 — 


water, and in substantial steady- -State reductions in downstream salinity with 


The effects of irrigation management on steady- -state ground- -water salinities 
_ have also been analyzed (6). As with surface waters, reduced leaching, at steady ? 
- state, may or may not reduce degradation of ground waters receiving irrigation 


_ 'Geochemist, United States Salinity Lab., United States ih of Agr., 4500 Glenwood — 
; Soil Scientist, United States Salinity Lab., United States Dept. of Agr., 4500 Glenwood 
= .—Discussion open until August st 1, 1981. To extend the closing date one month 
: written request must be filed with the Manager of Technical and Professional Publications 
4 ASCE. Manuscript was submitted for review for possible publication on January 17, : 
1980. This paper is part of the Journal of the Irrigation and Drainage Division, Proceedings — 
of the American Society of Civil Engineers, ©ASCE, Vol. 107, No. IR1, March, ee ‘4 
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drainage waters, depending on water type and hydrologic conditions. These 
_ evaluations were made for situations in which it was assumed that no salts, 
other than those present in the irrigation water or those derived from the dissolution 
of calcium carbonates or silicates in the soil profile, contributed to ground-water 
salinity. If additional readily soluble salts are present in the soil, or saline waters 
Dt. present in the drainage water flow path (as is often the case), then reduced 
_ leaching will always decrease the rate of ground-water salination. en Signy oN 
These previous studies have only considered steady-state conditions without 
‘considering travel times in either the unsaturated or saturated zones. Reduced | 
teaching substantially increases solute residence times in and below the rootzone | 
: _ (4) and thus influences the rate of groundwater salination. The purpose of this 
_ study is to demonstrate the important transitory effects of changes in irrigation » 
management on ground-water salinity, including consideration of solute transport "s 
in the unsaturated zone as well as chemical precipitation. 
_ The effect of irrigation management on water quality can best be determined — 
by computer simulation. Although the chemical and physical components of 
the model have been independently tested and verified, at the present time 
is not possible to verify the simulation with data from a real -ground- -water 
basin. Accurate historical data on leaching fractions, irrigation, and groundwater 
quality and soil CO, levels are not available. It is our belief that a major reason 
= quality implications of irrigation management have been neglected is that. 
_ its effects are observable only over long time periods (often decades or longer). 
Because of the considerable lag between implementation and observable results, 


“any changes i in water quality could also be caused by changes in in other var riables 
within that time frame. 


“time, and irrigation. water type « on n ground-water quality will be demonstrated 
; with a hydrologically very simple transport model. Although the approach taken 

here does not apply to any particular ground-water basin, the calculations serve 

to illustrate the underlying physico-chemical principles. Changes in soil ss 


soil layering, and depth to water table may therefore affect the numerical results 


of the analysis, | yet they \ will not alter the qualitative - conclusions | to be draw: 


Figure 1 gives a schematic diagram of the hypothetical ground-v water basin 


used in the analysis. _A water table is assumed present at a depth of 20 m 
below the rootzone, while the unconfined aquifer is | assumed to be bounded 
below by an impermeable layer at a depth of 50 m. . The ground-water basin 
4 managed to maintain the water table at a depth of 20 m, with no mass 
flow into or out of the basin. Because of evapotranspiration, this requires importing 7 
surface water into the basin at a rate equal to the net evapotranspiration rate, 
aaa minus rainfall; E is taken as 1.00 m/ yr). Water for leaching ve 


is pumped from the underlying groundwater through a series of wells which 7 


‘ drain water uniformly from the 30 m thick saturated zone. Thus, the rate a 
| pumping per unit surface area, Q, is equal to the average drainage flux, D. ‘+ 
_ ‘The spatially ‘distributed wells in the basin are assumed to be sufficiently close _ 


~ together so that, on the average, horizontal flow is not significant in the saturated 
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“GROUND. WATER QUALITY 
one. The basin is also assumed to be irrigated uniformly. Then no important 


areal ‘effects are present, and the ground-water basin can be treated asa 7 
one- ~dimensional vertical apn. Such a hydrologically simple system lends itself 


7 simulation of an actual gound- water basin may require the | use of a more venti: ; 


4 two- or three-dimensional model, this will depend o on the 


the vertical flow component, which i is ; examined i in ‘the study sodarsels ‘hao 
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Representation of Simulated One-Dimensional Ground-Water | 

Basin: = Net Rate (Set at 1.00 m/yr), S = Rate at Which 
= Drainage Rate below Rootzone (0.667 m/yr 

for L = 0.4 m/yr and m/ye for L L= = © ond of ‘from 

Two different leaching fractions (L) are considered: L=0.4 and L = 0. l 


| management ¢ optimizes water ‘efficiency (0. 1). The leaching | efined 
od “as that fraction of the applied water that leaches out of the ar: can 


_| 
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in . which Ss represents the rate at which water is imported per unit surface — 
area (S = E). In the present case, D = = Q : and S = 1.00 m’ >/m? yr. _One 
may calculate from Eq. | that the drainage flux for L: = 0.4 equals 0.667 m/yr, re 
and for L = 0.1 equals 0.111 m/yr. Total irrigation rates are then 1.667 m/yr 
for 0.4 leaching and 1.111 m/yrforO.l leaching, 
Water and salt movement are simulated with a slightly modified version of 
4 the one-dimensional single-ion saturated-unsaturated transport model described 


in detail elsewhere (8). The basic transport equations however will be given ” 


here. This study analyzes the effects of different irrigation management strategies _ 
on resulting vertical salt distributions over time periods of several decades. — 


: of drainage water over long periods of time, the use of a a steady state rather 4 
a transient is justified. The appropriate steady-state flow 


in which h = pressure head; K = hydraulic conductivity; and x = depth below 
the rootzone. The sink term Q(x) accounts for the withdrawal of water from — 


the saturated zone, and is given by: 


zero in the zone, then remains constant between the re 
table and the bottom of the aquifer. The following relationship between the 


—_— moisture content (8) a and the pressure head (h) was used for the | 


in which 6, (= 0.15 m >/m? ) = the residual _— water content; t; 6, (= 0.40 m 40 m >/m? ) 
= the saturated water content; a = —0.01; = 2; and m = | — I/n. The 


"parameter values given here are re typical for a a fine sandy soil. The howe 


ich K.: = 0.208 mm/s (or 15 is the sat hydraulic 
and 6 = the dimensionless moisture content: 
“Figure 2 gives a graphical representation of Eqs. 4 and 5. It should be kept 
in mind that the selection of the soil-hydraulic properties will have very little © 
effect on the results since the volumetric flux is fixed d by the rate of see . 


7 
7 / This makes it possible to consider only steady-state water flow. Earlier studies 7 
Wiereng and ui nd Reeves (1) suggest that to dict th lalitv 
. (2) 
a 


 egtene ad pumping rate. . These rates are in turn fixed by the aquifer depth 
_ and leaching fraction. Equation 2 is solved, subject to the boundary conditions - 


in which X (= 50 m) = the depth to the impermeable layer. Equation 7b represents 
a no-flow boundary condition. Note that the boundary flux at x = 0 is a 


+1000, 


Die 
levies 


200 28 30 


MOISTURE CONTENT, 6 fem$ fem?) 
FIG. 2. Relationship between Pressure Head (h), Volumetric Moisture Content | (8) | 


and Hydraulic Conductivity (K ) for Conditions Simulated us 


s the average 
- flux (D). The boundary at x=0 corresponds to the bottom of the rootzone; / 
the thickness of the rootzone itself is assumed to be be negligible compared v with vith = 


in which c = the ire salt concentration; Z = the dispersion coefficient; and 


7 
4 
| 
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q 
in which a = 0.0116 mm 2 A= 000 n mm; and where Zand q are expressed 
= in mm’/s and mm/s, respectively. The parameter \ is often referred to as_ 
_ dispersivity. The 1,000-mm value chosen for \ is much higher than is generally 
the case (10 mm-30 mm) in most soil physics studies. The high value of \ 
is necessary because of the Steady- “state ate approximation of the water flow equation. - 
This approximation results in an average downward flow of water. of only a 
few millimeters per day. Because of temporal variations in the irrigation and — 
evapotranspiration rates, both during the day and between days and seasons, 
the actual flow velocities of the water will be much higher and often change 
directions. These flux. variations in time will result in an increased — 


apparent dispersivity of the soil. “Additional dispersion effects are caused 
_ the fact that the basin is treated as a one-dimensional vertical system. Actual 3 
evapotranspiration and irrigation rates may not be uniform over the basin; 
_ irrigation frequency is often also not uniform while some portions of the basin 
- may not be irrigated at all. Also, the wells may not be distributed uniformly 
_ over the basin. These nonuniformities undoubtedly will lead to some lateral 
flow components, especially in the saturated zone, thereby increasing the apparent 
_ dispersion of the salt along the vertical dimension. The salt concentrations, 
therefore, must be viewed as areal averaged quantities. - — ae 
; 4 _— and boundary conditions imposed on Eq. 8 are: 


ac 


+ 


i. which C, 
_ The otivaiines material is taken to be a fine iil soil with no ion exchange 
potential. Two different water types are considered in this study. In each case 
® initial salt concentration in the soil profile (C,) is taken to be the same 
-s the concentration of the imported water. In the first example, «od is taken 


to be 8.52 meq/L with a composition “resembling that of Colorado River Or 


Si Grande River water equilibrated with CaCO, at 0.03 atm CO, pressure. 


Thus, this water is initially saturated with CaCO, (‘‘Type 2” water). Its 
composition it in per mer is as } Ca, 6.87; Mg, 0. 60; Na, 


The dispersion coefficient Din Eq. 8isevaluated by 
| 
7 
4 7 
—(X,t)= 
Cemical equilibrium calculations, based on a procedure cimiles to thet of | 


ND-WATER ALITY 


Ref. 4, showed that the salinity of the drainage water (C,) can _— ay 
related t to > the Salinity | of the irrigation water C, ) given | the foregoing. This 


composition ‘selected here. In addition, when Ca is not equal to HCO, (in 
milliequivalents per liter), no generalized relationship can be proposed between 
the concentration of the irrigation water and the concentration of the drainage — 
water. This is because waters of the same salinity but with different Ca 
concentrations can be derived by mixing various S proportions of different drainage — : 
waters generated from the same original water. In that ¢ case, the flow equation 
and the chemical equilibria model would need to be combined into a multi-ion : 
model and solved for each of the ions in the model. Again, this was not 

- the case for the water chosen in this paper, since Ca = HCO,. oe “a ome ody 

For 0.40 leaching, the linear relationship is 


For 0.10 leaching, | the relationship i is 


tation and concentration, it can be neglected fc for our analysis. 


For 0.40 leachin we hav e 


C,=0.1C, ,+0.9C, 0.1) . 
in which C, (= 8.52 meq/L)-= the concentration 1 of the imported water, a ; 


Co= the average concentration of the pumped groundwater ied ae 


LC into Eq. 12 allows Cat to 0 be expressed in terms 
Cy varies linearly with C,. 
analysis assumes that no chemical reactions occur below the bottom _ 
of the rootzone. The partial pressure of CO, in an aquifer below an irrigated 
- area with no recharge other than drainage water will be approximately the same 
as tha that of the bottom of the rootzone, if no CaCO, precipitation or dissolution — 
occur occurs when mixing waters of different composition. Initially, calcite-saturated , 
_ solutions can be supersaturated or undersaturated when mixed. This effect is 
explained in detail by Wigley and Plummer (12). As shown in their analysis, 
for fixed Poo, and temperature, this effect is largest when chemically dissimilar | 
waters are mixed, €.g., water where Ca > HCO, mixed with water where | 
Ca. < HCO,, and results in supersaturation. This effect is not important in 
our simulation. The other important potential cause of nonequilibrium is due 
to the nonlinearity of activity coefficients with ionic strength. This nonlinear ~ 


i 
— 
| 
i 
i CUUALIVU y at tiie ol CaCO, precipitated Is a linear function 
a _ of the concentration, with constant chemical precipitation for each leaching — 
| 
if im 
| 
» | P 


causes an in mixtures with different ionic ‘ 


__ strengths. The effects, however, are not substantial in the present — 


a because subsurface mixing occurs only with waters that differ at the most one 
_ order of magnitude in concentration, not between waters with orders of tinged 
differences in concentration. We assumed, therefore, that no substantial chemical a 
precipitation or dissolution reactions take place below the rootzone in the — 
- simulation described here. The ionic strength effect is more important when 
h. groundwater is mixed directly with dilute surface water; however, this is | 
accounted for by precipitation in the rootzone (Eq. 12), 
-... other water type considered is that of a water containing only Caso, 
_ The ground-water concentration at t = 0 is taken now as C, = 5 pay 
the same as the concentration of the imported water at all times. yor from a 


chemical equilibrium calculations, it was found the concentration of the 
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by 
FIG. 3.—Moisture Content (6) and Volumetric with for 0.1 0. 4 


> ple 
| 
+ which C,, as before, is given by Eq. 13 and Eq. 16 for 0.40 and 0.10 le 4 7 


¥ 

The value of C, ‘cannot exceed 31. 1s gypsum 
will precipitate at 31. 15 meq/L in a pure CaSO, system (25° C). Substituting — 
; Eg. 15 into Eq. 176 shows that the limiting concentration occurs immediately 
an the start of the simulation experiment, when the leaching fraction | equals: - 


bes:  Steady-s state distributions of the soil water content and the volumetric | a 


for both fractions are shown in Fig. 3. The volumetric: flux lux in the 


: 


oO 


unsaturated zone equals the rate at which water is pear out of the saturated 
7 zone: 0.667 m/yr for 0.40 leaching, and one-sixth of this amount or 0.111 m/yr _ 
. for 0.10 leaching. As shown in Fig. 3, -soil-water content in the ay 
zone is 0.238 for 0.40 leaching and 0.210 for 0.10 leaching for the given fluxes 
a The volumetric flux is constant in the unsaturated zone and decreases linearly 
from the water table to the bottom of the aquifer where it equals zero; this 
\ ll results from no water being pumped from the unsaturated zone and from uniform _ 


The mean residence (or travel) time of the salt in zone, 


+ which 20 = the thickness of the unsaturated zone. From Eq. 
18 it follows that the mean residence time in the unsaturated zone is about a 
a for 0.40 leaching, and 38 yr for 0.10 leaching. Thus, it will take more ar 


han five times a 2 for the salt to travel from the rootzone (assumed to 


CONCENTRATION, 


tion Versus Depth) Obtained after irrigating for. 
100 yr with Type 2 Water (CaCO, oo 
8.52 meq/L and Leaching Fraction of idl 


be of negligible thickness) to the water table for 0.10 leaching, as prone 
Figure 4 shows salinity profiles obtained after | irrigating for 2 yr, 4 yr, 10° 
yr, 20 yr, 50 yr, and 100 yr with Type-2 (Colorado River) water and 0.40 leaching. — 
The concentration distributions are quite diffuse at all times, partly because 
q of the relatively high value of the dispersivity used in the calculations, and 
partly because of the initially small concentration difference between groundwater 


and drainage water. This concentration difference will increase slowly after 
5 yt, when the « concentration front reaches the This occurs 


| 
Vy 
5.—Salinity Profiles (Concentra 
s 4 2 yr, 4 yr, 10 yr, 20 yr, 50 yr, and Gg 


GROUND- WATER QUALITY 
- a ow years earlier than the mean residence time (ry yr) because of the ‘effects bl 


= dispersion on the concentration front. At that time the average concentrations | 
in the ground water (C,) begin to increase, leading to an equivalent increase _ 
a in the concentration of the irrigation and drainage water (Eq. 13 and Eq. 12, 
respectively). The concentration distributions, however, remain quite uniform 7 
versus depth, especially in the unsaturated zone (Fig. 4). Note that the concentra- 
tion front reaches the bottom of the aquifer after about 50 yr.” 
_ Similar distributions for 0.10 leaching are shown in Fig. 5. While the concentra- 
tion i in the unsaturated zone now increases much faster than for 0.40 nadie, ) 
the concentration front itself moves much slower towards the water table. The n 
= edge of | the front reaches the ground-water table after about 20 yr, 
realy 17 yr earlier than predicted with Eq. 18 if no dispersion were em. 
: ‘The most striking differences between Fig. 4 and Fig . 5 are the gs and [ 
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.4 and 0.1 Leaching: Initial Ground-Water ‘Concentration, C, is 8. 52 52 meq/L 


, 4 locations of the solute | fronts. Low leaching causes much more salt to be : stored 


in the upper parts of the soil profile, leading to a slower salination of the 
ground-water system. Especially the lower part of the saturated zone 
remains relatively free of salt for a much longer period of time. | 
Since the amount of "imported water ‘is just equal to the net evapotranspiration 
’ ed with CaCO,, the total 
amounts of salt precipitated in the 1 rootzone must be the same for both leaching ~ 
- fractions. This can be seen by comparing Eq. 12 and Eq. 13. Precipitation — 
accounts for a loss of 61.8 meq/L for 0.10 leaching and 10.3 meq/L for 0. 40 j 
5 leaching. After correcting for differences in drainage volumes (six times higher — 
q for 0 0.40 leaching), we obtained the same mass of salt precipitation in are ® 


iy case. Therefore, the total amounts of -_ stored in the soil profile, including - 


| 


for 0.10 leaching. Figure 6 shows that this is indeed the case. The difference 
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both unsaturated and saturated zones, must also be the same. This \ will essentially a 
also be true for the average concentration in the soil profile, with small differences 
_ caused only by the different soil-water contents in the unsaturated zone and 
thus slightly different total volumes of water for the two leaching situations. — 
_ The most important difference between high and low leaching is a more uneven 7 


distribution of salt for low leaching. Because less salt is stored in the unsaturated 


zone, the average ground-water salinity, C,, should be higher for 0.40 leaching 


in ground- -water concentration for the CaCO, saturated (Type 2) water ~ 


Lia FIG. 7.—Salinity Profiles (Concen 
2 yr, 4 yr, 10 yr, 20 yr, 50 yr, a “100 yr with Type 3 Water (C280, 
Initial Concentration of 5.0 and Leaching Fraction of 0.4 
is only about 2.5 meq/L, and remains f airly constant in time. This concentration 
_ difference, however, is expected to increase as depth to the water table increases 
because the meen residence time of the salt in the unsaturated zone 1 increases — 
Calculated salt distributions for Type 3 waters (those capable of precipitating — 
gypsum) are shown in Fig. 7 and Fig. 8 for 0.40 and 0.10 leaching, respectively. 
q The curves for 0.40 leaching (Fig. 7) are shaped similar to those for 0.40 leaching 2 
Bi3 ‘with CaCO, saturated water (Fig. 4). This was expected because the travel 
times in the unsaturated zone are exactly the same for be both water types. The 


3 

— 
- jive @ 
| 


eine for most time intervals in Fig. 7. This i is due to the greater "amount 
of precipitation for the gypsum precipitating water (Fig. 7? than for the CaCO, © 
precipitating water (Fig. 4). The concentration of a pure gypsum water, further-_ 
more, cannot exceed the concentration of a saturated pure gypsum ‘solution 
— GI1.15 meq/L). This upper limit on the salt concentration in the profile leads - 
to a constant concentration between the 0 m-37 m depth after 100 yr (Fig. 2 
7). Gypsum saturation of the drainage water was reached after 75 yr of irrigation a 
for 0.40 leaching. This point was reached immediately for 0.10 leaching, leading , 


we 


2 yr, 4 yr, 10 yr, 20 yr, 50 yr, and 100 yr with Type 3 Water (CaSO, ~_ with 
initial Solute ute Concentention of 5 5. and Fraction of 0. 


to considerably more gypsum precipitation in the rootzone. concentration 


bi 
ig ‘8. —Salinity Profiles (Concentration Versus Depth) Obtained after irrigating for 


because of the diffusion-~ -dispersion. effects o on the salt distribution (Fig. 
Increased precipitation from high to low leaching causes the differences in cae: ey f 
distributions between the two leaching fractions to become more pronounced 
for Type 3 (gypsum precipitating) water. This is more clearly shown in Fig. 

9, where the omnes ground- water salt concentrations are plotted v versus } time. 


VAUTY 
‘ 
| 
| 
| 
=—linearly from 5 moq/L (the initial concentration) after shout 6 yr, to nearly 


4 
4 
4 


27 meq/L after 80 yr. The increase in concentration a deat down after 75 yr, 
since it can never exceed the gypsum saturation value of 31.15 meq/L. In i 
_ comneaat, the average ground-water concentration for 0.10 leaching does not 

appreciably increase during the first 30 yr of irrigation. After 30 yr, the rate 

of increase in ground-water salinity is also much slower than that for 0.40 
leaching. The average concentration reached 25 after 75 yr for 0.40 
leaching (Fig. 9), whereas this point was reached only after 200 yr for 0.10 ‘ 

leaching (not shown in figure). The maximum difference in ground-water salinity 

_ between high and low leaching was reached after 90 yr (28 meq/L versus - 

both cases will in saturated ground water; 


TYPE 3) Water 


CONCENTRATION, Co (m 


30 


FIG. 9 ee ae Water Concentration, & with Tene for Type 3 Water, an and = 
: for 0.4 and 0.1 Leaching Fractions: Initial Ground- Water Concentration, C,, is 5.0 


en erepeen state wee obviously will be reached much later for the lower leaching 
Any of the t type | given in this pe paper. is highly s site specific and depends 
upon the unique conditions prevalent in each ground-water basin considered. 
For the closed, shallow ground-water system considered here, it is apparent a 
- the benefit of low versus high leaching depends very much upon the type 
a water used for irrigation. Reduced leaching with Type 2 (CaCO, saturated) — Bs 
water has only a small effect on the calculated average ground-water salinity. 
- effect of low leaching would have been more pronounced if the depth 
to the ground-water table were greater than the 20 m assumed in the present 
calculations. The degradation | of the ground water would be delayed in time — 


Bi 
\ 
| 
| 


= 
and, additionally, a la a larger quantity o of Gasolved salts \ would then be | be stored in in 
“the unsaturated, rather than the saturated, zone. The main effect of increased — 
_ irrigation efficiency in the case of Type 2 water is the increased travel time 
in the unsaturated zone. This can be of considerable benefit whenever the depth 
to water table is large or if the irrigation water salinity is greater than — 
ground-water salinity. The considerably slower travel times for low leaching 
means that the saline drainage front takes much longer to mix with and degrade — 
the higher quality ground water, 
From Figs. 7, 8, and 9 it is evident that low leaching can be of considerable ‘ 
- benefit in reducing salinity when Type 3 waters are used for irrigation. The 7 
degree to which the benefit occurs is proportional to the amount of gypsum 
precipitated in and below the rootzone, and, therefore, depends wpon the 
composition of the irrigation water. The greater the percentage of precipitable a 
Salts in the applied water, the greater this benefit will be. For example, .= ai 
4 predominantly NaCl | type water with lesser amounts of Ca and SO, will not 
| ¥ result in much gypsum precipitation, and, thus, not reduce salinity to any great ‘ 
Pa: extent. The present example, that of a pure CaSO, system, represents a case 
4 of maximum benefit. The main reason for this is the increased precipitation 
with reduced leaching. As the proportion of soluble salts (or imbalance in Ca, 
HCO,, or Ca, SO,, or both) in the irrigation water incre increase, the benefits| of 
an improved irrigation efficiency will decrease. cnt 
_ In any actual analysis, the effects of ion exchange should also be taken into 
account. Ion exchange decreases the rate at which the salt front of the adsorbing — 
ion moves through the soil profile, and also may result in the precipitation: 
of additional calcium carbonate and gypsum. - Unfortunately, the effects of i ion 


capacity wal soil. For example, the exchange of Na or Mg for Ca does 
not necessarily induce precipitation. One may expect that the quantity of 

: exchangeable Ca precipitated will be larger for reduced leaching when HCO, 
>> Ca or SO, >> Ca. If, on the other hand , Ca >> HCO, and Ca >> 
 $0,, then relatively less of the exchanged calcium will precipitate. . In a recent 
study, Jury et al. (2) showed that the exchange of Na and Mg for Ca doubled 
the predicted quantity of salts precipitated in the top 150 cm of a particular = 
_ soil profile that was initially high in exchangeable Ca. They also showed the 7 
importance of the water-uptake distribution in the rootzone on the time necessary : 


to achieve steady-state salt flow in the soil rootzone. These latter effects are 
P robably less important for the type of analysis given in this paper, since the 
rootzone generally comprises only a small portion of the unsaturated zone. _ 
4g fact, in the present analysis, the spatial effect of the rootzone was omitted > 
altogether. If the rootzone were included in the simulation model, it would» 
serve only to increase the travel time in the unsaturated zone and would not 
significantly affect the calculated salinity of the ground water. ka 
: This paper considers a ground-water basin which is a closed hydrological 
unit, and in which the ground-water table remained at its initial position. If 
the hydrological system were to be altered by not introducing water, or by » 
significantly reducing the amount of imported water from outside the basin, — 
then overdrafting would be simulated, provided of course that the entire basin 
is irrigated with the same amount of water. . Then the salinity front would never 


| 
| 


__ Many basins are underlain by saline ground waters, with better quality \ waters 


— 
the re rate of propagation 1 of the solute front through the unsaturated zone. me. in 
this instance, salinity buildup in the ground water is not a factor. Maximum _ 

utilization of the available groundwater, however, would still require that high 


efficiency irrigation practices be applied. 
In many cases, the ground- -water basin will not be hydrologically — 


“surface: or subsurface flow system. The analysis will then become. more 
complex and probably not lend itself to a one-dimensional analysis of the type 
given in this paper. Nevertheless, we can still make some qualitative remarks, 
partly based upon the foregoing analysis. For example, if the applied irrigation _ 
water is imported surface water, then high leaching will probably result in a oe ; 
lower ground-water salinity because a greater volume of water and, thus, salts 

are being discharged from the basin ot oc) 


‘being located at much greater depths. The objective in this case would be to 
_ cause as little displacement of the saline water as possible. Mixing of the ‘saline — 
ground water with the ‘better quality ground water could result in a rapid 
degradation of a valuable water resource. The important variable now is not . 
the drainage water composition but rather its volume. Low leaching here would 
Se An optimal irrigation management strategy also requires judgment as to which 
_ water resource is the most valuable. For example, in regions with inadequate 


_ surface water storage, but with a potential for good quality ground water, priority 
should be given to maintaining or improving the ground-water resources. If, : 
however, the ground water is unsuitable or only marginally suitable for use, 
then the downstream river salinity would probably by the most important 
consideration in in deciding which irrigation practice tc to » adopt. ia 
A one-dimensional numerical analysis of water and salt movement in a 50 _ 
deep unsaturated-saturated soil profile showed the of leaching 
= and water type on ground-water salinity. . For CaCO, saturated water, 
the effects of reduced leaching are relatively minor, expecialy if ‘the depth ee 
to the water table is small. With i increasing thickness of the unsaturated zone, 
the mean residence time of the salt in the unsaturated zone increases, leading -_ 
_ to more salt storage in the upper part of the soil profile and less in the ground = 
water. The average ground-water concentration in that case will be smaller — 
For waters approaching gypsum saturation, improved leaching efficiency will 7 
significantly improve ground-water salinity, even in shallow aquifer systems. — 
For the case examined herein, salinity levels were up to 14 meq/L lower for 
0.10 leaching as compared with those for 0.40 leaching. The long-term benefits _ 
of reduced leaching are the increased storage of precipitated gypsum and soluble 
salts in the unsaturated zone, and a lower salinity of the underlying groundwater ; 
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‘The following are used in 


“concentration of drainage water leaching rootzone, in milliequivalents 


initial groundwater. concentration, in milliequivalents per liter; 


average concentration of pumped groundwater, in milliequivalents per 


dispersion coefficient, in square | millimeters per second; 


year; 


= pressure head; 


saturated hydraulic conductivity, in millimeters er r second; 


Dall Burdens Kesulling Irom igauon With Kiver Walters in ine Weslemm United 
sou characteristic parameter used in , IN Square m1 imeters per 
= 
a 
4 
q Co = 


= 


= 


= : fraction of water ‘that leaches out of 


m = soil characteristics parameter used in Eq. 4; 


soil characteristics parameter used in Eq. 4; 


= volumetric flux, in millimeters per second; 
Q = rate of pumping per unit surface area, in cubic meters per square meter | 


: S = rate at which water is imported into basin per unit unit ‘Surface area, in 


depth of impermeable la cr, ¢ 
thickness of saturated zone; 

soil characteristic | parameter used in Eq. 4 
8, residual moisture content; 
© = dimensionless moisture content; 
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Economic MopELInG 


By Henry oO. and Robert W. Hill,? A. M. ASCE 


 Thei increasing use of water for i irrigation, industry, and municipal water supplies 
_ has led to an increasing need for careful planning for the usable water resources. Ps 
Efficiency in water use requires skilled planning and careful management of 
water. Efficient use considerations in a dynamic system such as a river basin 
implies that the physical and ¢ economic systems | be described with sufficient 
accuracy to quantitatively and qualitatively predict the system- -wide effects of _ 
_ depletions resulting from water use anywhere in the system. Unfortunately, | 
~ such depletions for irrigation and industry concentrates and adds nondegradable 
_ substances that produce a decrease of water quality. Thus, s every hydrologic 
system, each upstream u use ¢ has some e effects o on n the ; quantity, » quality, and timing 


: benefits of var various - upstream 1 changes or ‘management alternatives can lead to 
increased use and better efficiency of available water resources within the basin. a 
An appropriate description of a water resource system, therefore, includes wad 

= hydrologic system, the salinity flo flow system, the economic system, and those — 
functions which relate them. | rez 
_ This paper presents a 
be combined with an existing “hydrologic salinity model in n order to bl 
‘Management effects of sequential water use in a river basin. The hydrologic, “y 

2 and economic systems are closely interrelated in any water resource 

‘Shee bere planning is difficult if the three systems are analyzed 


= 


 Lect., Dept. of Agricultural Engrg., Univ. of Ife, Ife, Nigeria. = = | 
~ Note.—Discussion open until August 1, 1981. To extend the iil date one month 
4 written request must be filed with the Manager of Technical and Professional Publications, _ 
; ASCE. Manuscript was submitted for review for possible publication on June 3, 1980. : 
_ This paper is part of the Journal of the Irrigation and Drainage Division, Proceedings 
of the American Society of Civil Engineers, OASCE, vel 107, No. IRI, March, 1981. 4 
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FROM OTHER UNIT ECONOMIC | 
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where 


“a 


wegen 


Salt Tolerance Function as Link Between Salinity- 


10 
| 
a 
q | 
; 


othe hydrologic, salinity, and economic flow systems are dynamic in nature 


_ and can be related through the concepts of continuity « of mass and m momentum. 
Fens of momentum may be considered negligible in this case since velocities - 
g low. A proper accounting of the physical hydrologic quantities is obtained 


TABLE 1 1 _—Salt Tolerance of Crops (7) 


yield decline 
A, decrease per 

percenti- § | salinity beyond 
meter (decisie- threshold B, Sait 
mens per meter) ) | asapercentage tolerance 


‘Sugar cane 
Tomato 
ed Wheat 


2 


as as water is translated or routs eddie a the sys system with respect to both space 
and time, described by the expression (adapted from Ref. 5): 


+ Qu + — ET + Q,, + Q,0) = 
in which P.= the precipitation on ‘the: area; Q,, = the total surface inflow; _ 


7 Q, ie = the subsurface inflow; ET = the evapotranspiration from the area; 0. 
4 the total surface outflow; Q,, = the subsurface outflow; and AS = the change — 
‘in water storage including snow, soil moisture, surface reservoirs, and ground 
a water. ‘With the exceptions of precipitation, snowmelt, and evapotranspiration, 
all of the water quantities described as the hydrologic model in Eq. 2 have 
quality parameters associated with them. If all these quantity and quality 
‘components making up Eq. 2 can be identified, then the ; general salt salt simulation — 
“model cz can be in form as: am 4 


2 «eer 


subbasins in a sequential manner moving downstream. The model providesa 
| 
Crop 
— 

| 
i 


= 


0... = => + + Cin 24 AS, C 

in which Ci = the salt concentration of outflow water; C= the salt concentration 


isj 


of surface ‘source, Js Cue = the salt concentration of ground-water source, k, 


in outflow; and C,., = the salt concentration associated with tana element, p. a 


ComeininG Saunity ano Systems 
4 
Careful consideration should be given t to the nature of the re iationships existing — 
among” ‘the hydrologic, salinity, and economic systems within a river basin. 


The hydrologic and catinity flow systems are physical and quantitative pee, 


MBINED ECONOMIC SIMULATION MOCEL 


SOW. MOISTURE DATA 
“WATER HOLDING CAPACITY 
INITIAL, SOIL MOISTURE 


SALINITY. para AND 


INPUT TIME PERIOD WEATHER _ PARAMETER: 
PERCENT | SNI, SNK, 3M, TPR com, CNA 


MOURS, ETC 


ok 


CALCULATE WATER 


DETERMINE MOISTURE 
STATUS DE! 


AND 


TE RELATIVE 


BY THE 


—Flow tow Diagram of Combined Economic Simulation | Model 
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while the economic model is simply a quantitiative phenomenon. When combined, — 


- these three systems form an integrated unit since there exists a strong interrela- — 
_ tionship between the physical and economic systems. Fig. 1 presents the basic 
_ components of the hydrologic, salinity, and economic systems and the relation- _ 
ships existing among them. The economic system includes the crops produced = 
on irrigated farm land and their relationship to the hydrosalinity system, but 


does not directly consider livestock Or municipal uses. The basic economic 
unit of the study isariversubbasin, 
_ The connecting links between the hydrologic, salinity, and economic flow 
systems of an agricultural enterprise are dependent upon such factors as water 
availability (supply), water 6 a (demand), ‘Salinity levels of f irrigation 


(KANO RIVER 
PROJECT 
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4.—Map of Hadejia River Basin 1 Showing Subbasins Used in Study boa 
water and soil solution, crop salt tolerance level , and production per unit of | 
that influence production i in agriculture. Production is also a function of r agg 
‘ment, capital, labor, crop variety, weather, and soil type and fertility. yo 
maintaining these other factors at relatively constant levels, production (yield) 
can be estimated at known water and salinity levels. The link between the 
: hydrologic and the economic systems is the production function for each crop 
A production function determines the relationship between the variation 
in yields of several crops resulting from a variable input of water as influenced 7 
by salinity levels of the soil, while all other factors are assumed held constant. — 
The link between the salinity and economic systems is the crop salt tolerance _ 
function later in this pat. a schematic diagram of 


water consumed. However, water and salinity levels are ‘not the only factors 


a 
— 
= 
ii 


wy the hydrologic, salinity, and economic a are linked together by their 

Production Function.—An estimate of the production fu unction is the re relationship § 
between the crop yield and the seasonal evapotranspiration of the crop, for 

_ the assumed constant or near constant levels of fertility, management, and other 

- conditions. The concept, as used by Doorenbos et al. (2), that relative yield — 
| decrease i is a function of relative evapotranspiration deficit through an ae 


= the crop yield: = the or potential 
m 
crop yield; k, = the crop yield response factor; ET, 
tion; and ET,, = the maximum or potential evapotranspiration. The values of | 


Y,, and k, are given for various crops in Ref. 2. The variables, ET, and ET, 


a all terms are previously defined. Eq. 4 is valid a ie water iis 
up to 50% or ET,/ET,, = 0.5. If water deficits exceed 50%, then the actua 
ie Y,, is assumed to decrease linearly from the point \ where ET /ET 
5 to the point where ET, /ET,, _=0. 0, and is expressed a 


in which Y,, = the actual harvested crop yield when water deficit exceeds 
50%; * = the actual crop yield when water deficit equals 50%; ET, 
— evapotranspiration when water deficit equals 50%; and ET, = ; 
evapotranspiration when water deficit exceeds 50%. Eqs. 4 and 5 thus become 
_ the production function for each crop and completely define ony yields at — 
all water deficits once the potential yield, Y,,, yield response factor, k,, actual 
and potential evapotranspiration, ‘ET, aan ET... Tespectively, are known for 


‘This approach i is good as initial input into the model and is especially vader : 
in countries where data on crop water use and yield relationships are rie 


sparse. Where complete and accurate data are available on crop water use, 
p 


_ Because of its simplicity, low di data requirements (only surface > air temperature. 
is needed), and applicability to most of the irrigated arid regions of the world, 
7 the modified Blaney-Criddle equation (9) has been adopted for the calculation . 


: of actual and potential evapotranspiration in this model. The equation is expressed — 


| 

| 
q 
| 
| 
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in which PET = = the potential evapotranspiration; k. = the monthly | crop oon > 


coefficient; and k, = the monthly temperature coefficient determined from the 
following equation = 


= the m daylight hours of the year, a asa percentage. 

TABLE 2.—Summary of Subbasins Seasonal Water Diverted and Salinity Effects 

on Net Returns, Water Years 1976-1978 (calibration) $= 


| 


Solution Salt Seasonal | a 
4 Concentration Net Returns, 


water di- 


2. In | In verted, in 
parts | millimhos | acre-feet 
i 
million | centimeter of land 


River 


Wudil- Dabi 


Dabi-Marke ke 
= 


Burum Gana River | 1976 | 468 oT) « 


P | 
1977 | 1,074, beg, 68 102 
Note: 1 mmho/cm = | ds/m. Multiply acre- -feet | per acre by 304.8 to obtain millimeters. 


¥ Gardner and Ehlig (4), among others, have gies that transpiration occurs 
at the full potential rate through approximately the first one-third of the available 
soil moisture range and that thereafter the 1¢ actual evapotranspiration rate lags 


- the potential rate. This phenomena is determined in the model by the following — ‘ 


> 


. aa 
subject 
| 
te 
Water | and 
4 4 | 
1976 | 629 | 098 | 126 1260 
&§ 1977 | 1.49 120, 
| 207 | 207 
&g 182 | 182 
which ET = the estimated actual evapotranspiration, in inches; CSM =the 
J 


MARCH 1981 


3.—Summary of Subbasins Seasonal Water Diverted and Salinity 4 


crop area | percentage ‘ a 
Nad = q ait 
| 1976 | 
4 
1978 
1976 
1977, 
1977 


1978 
 Hadejia 
1976 
1977 


Nigerian Naira = 1.77 United States dollars. 
Note: Divide acres by 2.47 to obtain hectares. Multiply acre-feet per acre by 3048 8 7 


_ critical soil moisture level below which ET is less than PET (in inches), determined 
_ during the calibration process of the model; SM = the soil moisture level, ; 
in inches; and SMC = the soil moisture capacity, in inches (root zone water 
storage available to the plant). The estimated actual evapotranspiration is then 
accumulated for the growing season and for each crop and is introduced as 
the basic input to the economic model fr: from the hydrologic system. 
a Crop Salt Tolerance Function.—Crop salt tolerance has usually been ‘expressed — 
as the yield decrease expected for a given level of soluble salts in the veel 
a= compared with yields under nonsaline conditions (1). The most common 
method of measuring salinity is to determine the electrical conductivity of 
ie extract, EC., , from the active root zone. The salt concentration of | = 
the soil solution was calculated | through | an accounting process at ss een toed 
model time period. This process can be expressed by: 
CSS = f(SSIR, SSWF, SSDP) ap 


i which CSS = the soil solution concentration; SSIR = the salt applied from 


- irrigation waters; SSWF = salt added from weathering or fertilizer application; — 
_ and SSDP = salt lost by deep percolation. Weathering rate of mineral is considered ’ ; 
significant only i in the ground-water system in contact with the oes — 


a 


| | 
133 | 085 | -36 | 
112 | 0.59 
-10 
A 


"Effects on Net Returns: Comparison of Calibration and Management {Option a 


Water andSalt 


percentage 


ad 


_ Maas and Hoffman (7) compiled and normalized all available salt tolerance 
data from an “extensive review of literature. Their findings were adapted into 
linear equations which relate the telative yield to soil saturation EC. 


= the relative yield for any given soil salinity exceeding the » thecal; ¥e 
as 3 a percentage; A = the salinity threshold, in millimhos per centimeter (in _ 
decisiemens per meter); B= yield decrease per unit salinity increase, as 
i percentage; and EC, = the soil saturation extract, in millimhos per centimeter — 
dn decisiemens per meter). The values of A and B were provided in the Maas 
j and Hoffman report | for several crops; some of their data relevant to this study 
were extracted and are presented as Table 
_ Egs. 12 and 13, subsequently referred to as s the crop crop salt tolerance function, — 
can be used in conjunction with the values of A and B from Table | to determine — 
the relative yield when the soil solution concentration is known. Thus, when | 
the yield of each crop is determined by the seasonal evapotranspiration, further 
modification (reduction) of yield is done the yield factor, if 


102 @ | 102 -17 
| 
| 


TABLE 4. .—Summary o of Subbasin Seasonal Water Diverted and Salinity Effects on 


per acre of irrigated - 


Change, | as a 
percentage year | Ml MNG percentage 
| 7) | 
1977 
Wudil-Dabi Dabi 4, 1976 
4 Dabi-Marke 23 1976 
1977 
1978 
1976 | 
1977 


 “MNGI * Results from i increas crop areas. 
[MNO II = Results from incre 
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: Combined Simulation Model.—The hy- 

: ttnale and salinity s systems have been programmed on the computer in Fortran 
Iv language by Hill et al. (5), and further modifications on these have been 7 

done by Huber et al. (6). The combined economic simulation model is presented 


asa conceptual mathematical model in Fig. 3 in which data o on crops—costs 


are combined into a a single working model. nai a tl 

_ The model, hereafter referred to as the Sequential Water Use Model (SWUM), 

: has as its main objective the evaluation of the marginal primary benefits of 
Pe by computing the incremental changes in net returns to the farm unit 
- the as a result of changes in the water supply. It also answers rs questions concerning 


_ ‘The model requires crop, soil, temperature, daylight percent, hydrologic and i 
salinity data, and calibration parameters as input (see Fig. 3). The model output — 

_ includes monthly and seasonal values of: potential and actual evapotranspiration, 
- irrigation water diverted and applied to each crop, and soil and — 
salt levels. It also computes the total irrigation water diverted and applied during 
§ growing season, the total water shortages and reservoir storage capacity 
"required to meet such shortages. Finally, the model calculates yields and estimated - 


net return per om area for each crop and for the entire subbasin. | 


| 
| 
Phreatophyte Area, in acres 
/ 
| 
-Burum 7,640 3,820 | 1976 | 0.73 
4 plus reduced phreatophyte areas. 
| 
? Note: Divide acres by 2.47 to obtain hectares. Multiply acre-feet per acre by 304.8 — : 
&§ 


Net Returns: Comparison o of Option | Option 
‘Cur 


_ |Change, asa; Change, as 


‘The relationships proposed in the p preceding : sections of were re applied 
Bs to the Hadejia River Basin, located in Northern ‘Nigeria. To account for spatial 
_ = _ differences in climate and basin characteristics, the drainage area was divided 
into five subbasins, as shown in Fig. 4. Data on streamflow, water quality, 
a meteorology, crops, and economics were taken or estimated from published 
¥ - sources. Water and salt outflows from the upstream subbasin became the inflow _ 
for the next one downstream in a sequential manner. ase ae 
_ The parameter set of the hydrologic and salinity processes must first be 
if determined by calibration before the sequential use model can be used. ng ae i 
a calibration parameters become inputs into the sequential use model. 
mathematical programming approach used to determine the parameter 
which describes the model involves selecting an objective function 


in which DIFF = the numerical difference between computed value and renga 
ug é value of a given hydrologic (streamflow) or salinity (e.g., tons of salt) d data. 4a 
Thus, the object of the calibration process is to adjust the various parameters” 
in such a to cause DIFF to approach this is 


&g 
| 

| | 
| 


Change _ Water Diverted, in acre- 
acre of irrigated land 


soil solution 

salt concentration 

Subbasins a percentage CALIB® 
0.96 

1.03 
1.37, 

| 133 
138 


*CALIB = Results from calibration runs. 
III = Results from drought simulation. moll 
Nigerian Naira = 1.77 United States dollars. 
a ‘Note: Divide acres by 2.47 to obtain hectares. Multiply acre-feet per acre ns 304. 4 
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the or is said to be calibrated. The model can then be used with relative ', 
~ confidence in subsequent management studies of the basin. its 
- The program’s capability to model water quantity and quality effects from # 
different management alternatives and subsequent predictions of net returns 
resulting from such | changes in ‘management is of tremendous importance to 
_ river basin planning and implementation. Some of the management options - 
+ Increasing crop acreages to projected values. 
2. Reducing phreatophyte acreages as a better management practice after” 


7 


_ 3, Simulated drought conditions. The main purpose of this trial was to observe 

the effect of low water supply and an induced highly saline soil condition on 

The results and analyses of calibration and management studies follow. A 
_ summary of all subbasins’ seasonal irrigation water diversions and consumption, 

and the net returns of irrigated land, is presented in Table 2. There is no general 
trend in the amounts of water diverted and consumptively used per unit area 


subbasin to subbasin; however the net returns per unit area show 


_ TABLE 5.—Summary of Subbasin Seasonal Water Diverted and Salinity Effects on 
| 
— 
hange, asa > 
092 
0.92 | 
1.01 
1.01 
if 
— 0.82 on 
— 
| 
| 
i 


Net Returns: Comparison | ot Calibration with Water Years 


pe 


Change,asa; | |Change, asa 


“MING | percentage | percentage 
9) (13) 


i 


ogy 
« 


; the se season with the highest | water diversions seems to have ‘the h highest c corres- 
: _ ponding net returns within each subbasin. The Hadejia Valley Subbasin with 
the highest diversions has the maximum net returns in all seasons. An examination 
os of the salinity effects in Table 2 shows that within the : same season 1 (water 
year), the concentration ‘tends to increase from subbasin to subbasin in the 
downstream direction except where tributaries with low salt concentration and a 
7 high discharge volumes dilute the river inflows. However, the predicted soil | - 
solution concentrations produced no significant yield 1 reductions and hardly any : 
decline i in net returns with the calibration runs. thee 
Increasing Crop ‘Areas.—With increases of some 59%- 100% (see Table 3) 
in irrigated land areas of the subbasins, the amount of water diverted per unit 
-*- showed reductions from the calibration run except on the most upstream a 
a Ss subbasin (Kano River). An evaluation of the net returns, with respect to water 
only, showed Kano River Subbasin with positive change (increases) i in net return 
a all of the three growing seasons. - Wudil- Dabi Subbasin had increases in net "i 


_ returns in two of the three growing seasons. The remaining subbasins had declines __ 
& net returns. This suggests that only the two most upstream subbasins could : 


provide adequate \ water for such increases in crop areas under the present _ 


_ Adding the effect of salinity imposes additional burdens on the water satay : 


situation. A further look at Table 3 shows no effect of salinity on net returns” 


= 
| 

| 

| | 
| 


in the However, from the second growing season 
. on the Hadejia Valley and Burum Gana Subbasins, salinity buildup resulted a 
ie — in declines in net returns per unit area. The reductions in net return were slight w 
in most cases; | - however it reached 11% | on the Burum Gana Subbasin in the 
Reducing Phreatophyte Aress.—The water areas were 
- reduced by 50% in all the subbasins as an expected reasonable improved - 
_ management practice after increasing the crop acreages to projected values. 
A summary of subbasin seasonal water supply and salinity effects on crop 
yields and net - returns for water years 1976- 1978 is presented in Table 4 for 
the management op options indicated. With 50% ‘reduction i in phreatophyte areas, 
the model predicted up to a 42% increase in water diverted per acre in one 
¢ of the subbasins; however some subbasins did not divert more water even with 
_the increased water availability. The logical explanation for the lack of increased 
__ diversion ‘might be | water r became available when not needed, i.e., more water ‘ 


to “surface and deep percolation losses. The downstream subbasins q 
responded more favorably to water availability by diverting more water. In 
all subbasins, increases in water diverted per unit area resulted in from 0%- 14% 
increases in net returns with the reduced phreatophyte areas (see Table 4). te 
The increases in net returns are more noticeable in the downstream autpecins, 
Sight or no increases at all are common in the upstream ‘subbasins. _ >. 

_ The additional effects of salinity were most noticeable in the downstream 
adteetii (Burum Gana River) where net returns were less for all growing seasons 
when salinity effects were imposed than with water alone. However, some = 

_ increases in net returns were observed when the phreatophyte areas were reduced 
after first increasing the crop areas. 

Simulated Drought.— With the of some growing seasons, the percent 
enm in water diverted going from calibration to the simulated drought condition 

- are negative (see Table 5), which, as expected, indicates less water availability 
for crops when drought occurs. On the other hand, the amount of water diverted - 
to Wudil-Dabi Subbasin in 1976 and 1977 growing seasons for drought simulation. 

were actually higher than for calibration. An error (under-estimation) in the 

‘measured canal diversion used for calibration may account for this discrepancy. 

In general, drought simulation resulted in up to 30% decline in water supply 

in some subbasins and more than 15% decline in most subbasins. As a —_ 

_ quence, there was more than 30% cut in net revenue for most subbasins (see 

5). ‘The downstream ‘Subbasins (from Dabi- “Marke to.  Burum Gana R 


¥ situation, the decline in net returns reached 43% on the Hadejia Valley Subbasin 
: in 1978. As was the case for water shortages, the downstream subbasins were 
the most limited when salinity effects were superimposed upon that of water, 


. 


mr. hydrologic and salinity model was combined with an economic simulation | 
_ model in in order to predict management effects of ‘Sequential water use and salinity 


| 
he addition of salinity effects aseravated the already had wate 


RIVER BASIN MODELING 
y subbasins in a river | The sequential use model thus developed uses 


"concentration to ) yields. In ‘essence, the outputs from the hydrologic (seasonal 7 
evapotranspiration) and salinity (calculated soil solution salt concentrations) 
models became inputs into the economic model. geet 
= For testing, the motel was applied to an actual hydrologic unit, the noe ar 


River Basin i in Nigeria. The practical utility of th the model was ‘demonstrated 


a "and salinity conditions. Management ru runs of the model indicate that for Hadejia 


1. Increasing the crop areas beyond those used for model calibration resulted 
4 in increased total net return for the entire Basin area; however net t returns ty 
(Per unit area were not as high as with the calibration runs. Declines in net 
_Teturns per unit area were due to water shortages and salinity buildup in some 
i * Reducing the areas occupied by phreatophyte, /c open water indicated that | 
water previously consumptively used by phreatophytes could be beneficially 


used d by the © Crops. However, such reductions i in phreatophyte. / open water areas” 


3. In the event that drought occurs, the model predicted water - shortages es 
-_ and salinity buildup» which eventually resulted in the decline of net benefit Ss 


totheentire basin, 


can n be useful in evaluating | the crop acreage e combination which will 
maximum overall benefit in evaluating project design reservoir capacities and — 
other facilities appropriate for the demand and available water. The model is 


also evaluating the relative efficiency of water with respect to 


sed 


confidence if the production functions and the salt tolerance functions had — 7 
defined more precisely. However, isolation of two factors (water and ay) 
as measurements of production in the complex process of crop growth and 
yield is very difficult and may not represent all the factors affecting crop growth. | - ) 
There are other factors assumed constant in the development of this model, 

which in reality are not. Such ‘“‘other factors’’ include soil fertility, infestation 
by pests and diseases, etc. Continued progress $ towards a more precise estimate 
of crop yields will require additional effort directed towards the “or wie 
of the effects of these these ‘other factors” on crop growth and yields. 
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B = yield decrease per unit salinity increase, as a percentage; 
salt concentration of ground-water source, k; 
Cy = salt concentration of surface source, j; 


salt concentration associated with storage element, p; 
concentration of outflow water; 
_ CSM = critical soil moisture level below which ET is less than PET; 
soil solution salt concentration (total dissolved solids); 
- DIFF = numerical difference between computed and measured values ues of 
given hydrologic or salinity data; 
EC, = electrical conductivity of soil saturation extract, in ‘millimhos | 
ET, = actual ET when water deficit equals 50%; 
= maximum or potential ET; 


total surface inflow; 


1 
| 
| | 
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soil moisture capacity y (root zone water er storage available t to plant); 


salt lost through deep percolation; __ On 
salt applied from irrigation waters; | 


added from weathering or fertilizer applications; 


= harvested crop yield; 


a = 
Y, = oma crop yield when water deficit equals 50%; 
maximum or potential crop and 
= = change in storage. 


thik 


‘SSIR = 

q | | 


+ 
AS 
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FOR SURFACE EFFECTS +t) 
By Tan D. Moore’ te: 
Surface sealing and ‘mulching can significantly influence infiltration and 
therefore surface runoff and soil erosion. These effects have been demonstrated mm 
by numerous researchers (8,9, 10,12, 13,14,18,21,24). Mannering (14), e.g., showed 
that in many instances infiltration rates on bare soil are only 20%-30% of the | 
. rates on protected soils. Surface sealing effects can greatly overshadow other — 
factors affecting infiltration on unprotected soils. Mulching and incorporated 


:* organic matter also have a considerable, , but opposite effect. These effects — es 
4 


i in infiltration modeling | is ; lacking. This paper shows how these “effects can 


be incorporated into the approximate two-stage infiltration equations proposed — 

based on Darcy’s law using a capillary-tube analogy. Swartzendruber (25) pointed — 

- out that the capillary-tube model was more restrictive than needed for Green 
and Ampt’s approach. The method assumed an initially uniform moisture content — 
ina ‘homogeneous soil. Philip (22) later derived the Green- -Ampt equation ci calling 

it the delta-function solution. Morel- Seytoux and Khanji (20) ) derived an equation ; | 

- the same functional form that considered air viscous effects of two-phase id 
ka but without many of the restrictions, including the capillary-tube ~ 6G nl 
oobi (4,5) used the Green-Ampt equation in his studies of infiltration “a 


So "Asst. Prof., Dept. of Agricultural Engrg., 106 Agricultural Engrg. Building., Univ. 2 
_ Note.—Discussion open until August 1, 1981. To extend the closing date one month, 
. written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on August 12, 1980. 
This paper is part of the Journal of the Irrigation and Drainage Division, Proceedings 
the American Society of Civil ©ASCE, ‘Vol. 107 IRI, March, 1981. 

ISSN 0044-7978 /81 /0001-0071/$01.00. 
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nonuniform soil profiles. The Green- -Ampt equation ras received considershle ; 

attention in recent years and although it is an approximate equation it has been 
_ Shown to have a theoretical basis, as well as measurable parameters. 

If the depth of ponding | is negligible, then | the Green- — } equation can be 


dF Ae GMA 


in which f = infiltrability; F= the infiltration = the initial 
x 
: moisture deficit; S, = the capillary drive at the wetting front; K = the hydraulic | 
conductivity in the wetted zone; and ¢ = time. This equation is a single- “stage 
— 
; infiltration equation because it assumes ponded conditions at all times. — —s 


A Ampt, developed an equation for | predicting infiltration volume | prior to surface 
ponding (assuming a uniform initial moisture profile and constant application 


‘in which F, = the: volume of water infiltrated at the instant of surface ponding; 
and J = the application rate (constant). Mein and Larson (16,17) then modified 
the Green-Ampt equation to predict infiltration after surface ponding (the second 

. ain of infiltration) by translating the the time sca scale. Thi This modified equation is = 


the time to infiltrate a volume ‘to F, under surface 
conditions (obtained by substituting into Eq. l and solving for Eq. 2 
can be solved explicitly for F, or t, and Eq. 3 | can be solved ‘explicitly { for ie 


the Mein- (GAML) Model. 


Derivation of GAML Mooirieo FOR SURFACE Errects am 


_ Modified Green-Ampt Equation.—Fig. | represents the conceptual soil profile 
= in the following derivation of the modified Green-Ampt_ equation. The 


| 
| 
| 
{ 
a ‘aah and characterized by a hydraulic conductivity, K,. Below this, 
i =the soil is homogeneous and semi-infinite with an initial moisture deficit, a0, 


If the depth of ponding is negligible, ‘then Darcy’s law be to 


a 
in which S$ = the capillary drive at the wetting front in the subsurface soil; 
7 bars (L + S,.) = the potential at the wetting front if the potential at the soil — 


| 


4 

«AIG. 1 —Conceptual Moisture Profile for Derivation of Modified GAML Equations 


— is is arbitrarily taken as zer time, 4 bas volume of water infiltrated 


| 


Tk, 


DEPTH BELOW SURFACE 


4 


Because f= dF/dt; 


| 
“> 
} 
a 
: 
|} 
to 
7 
| { | 
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nging Eq. 7 results in 
— + 46,5,,); and B= 


Merefor 


OF, on integrating» KK, In(A + K,K,F)=t+ c.. 
in which C=a “constant t of integration. The limits of integration a are j 

4 atr= ty and F = F at t= t. Applying these limits to Eq. ' 


(= ) 

(A8, 1) + (a8, +F 


) + (46, S..+F 


F, 


the 
. storage volume in the surface layer. "This equation was ; first presented by van 
- Duin (27) without details « the derivation, and later derived ved by / Asseed and and = 


of 
When L = L., , the Green-Ampt equation can be written as 


in which S,, = ‘the ‘capillary drive at at the wetting front in the surface layer 
of soil. The term fy in Eq. 10b can be evaluated by substituting F, for - 
in Eq. 10c. If a Toms then Eq. 106 reduces to Eq. 10c, written for uniform 
_ When L > Ly, the siete is computed by Eq. 3 which can be ‘rearranged 


and expressed as follows 


pre 
E+F-F, 
When Ll = 


(9b) 
add 
| 
: 
in which E = L,A0,(K,/K,); H = A0,(L, + S,,); F, = L,40,; andt,=time =f 
4 4 | 
(11b) 
, the 


NE EQUATIONS 


Ney 


A 


le“ 


soil profile at the moment of shown in Fie, 
Therefore, from Darcy’s law: 


top layer a at the t time ad surface ponding. From Eq. 1 


and it therefore follows thet 
Ry x 


= L,(A6, — + L, Ae, 303 

which results i in L,= (14c) 

Equating Eqs. ‘136 (Darcy’s equation) and 14c (continuity equation) yields 


— 
“un 
— 17 | 
= her 
4 
From Fig. 2we see that (14a) 


4 
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A6, 


L,A8, 


in ins E, F, and H, have the same definitions as in a J Modified Green-Ampt | 

equation (Eq. 105). Eq. 155 applies when L, > (i.e., when exceeds the 


storage volume in surface layer). When L, <L, the Mein-Larson equation 


by be written 
Modified Green- -Ampt Equation for Two-Stage Infiltration. —At time t, 
_F, and at time ¢,, F = F,. Therefore, for two- teenie eel when L > 


‘ 
L, and L, > L, Eq. 9a can be writtenas 


(16) 


“field saturated surface conditions (t’ is by substituting F, ‘for F 

- Setting F, to zero in Eq. 106 and solving for #). When L = L, and L, = 

* the modified Green- -Ampt equation for for two- vo-stags inf infiltration can written 


in which t'= obtained by substituting F, for F in Eq. 10c ae solving for — 
t. When L > L, and L, > L,, Eq. 176 cnn be used for predicting F from t = ty 
to t= 1,, and Eq. 10b used thereafter. If L, = 0, then 1 Eq. 17a reduces to Eq. 
Eqs. 17a and 176 are used to calculate the infiltration volumes after surface a 
4 ponding. The corresponding infiltrabilities are given by 
Parameter Evaluation for Modified GAML Equations. —Th —The modified GAML- 
equations cs can be be applied to cases in which the the surface layer is either: (1) More 


| 
| 
B+FK, | 
= 
| 


_ for the model (specifically S,, and S.2) requires a knowledge of the relationships 
_ between hydraulic conductivity and capillary suction (K versus S), and capillary 
suction and moisture content (S versus @). Parameters derived from these 
relationships have physical significance. An alternative, though notas theoretically _ 
Satisfying, i is to derive fitted infiltration parameters. Asseed a and 
(3) used this method to experimentally evaluate Eq. 17a. 
Ds _ When the ‘surface layer i is more pervious than the subsoil, the subsoil controls — 
infiltration rates. In such cases the hydraulic conductivities, K, and K,, for 
_ input to the equations should be either the saturated or rewet hydraulic conduc- — 
tivities in the two zones, depending on whether or not air entrapment effects 
are deemed important. The capillary drive terms, S,, and S,,, can be determined _ 
using one of several procedures (1,15,19). “Aggelides and Youngs (1) found that 


pervious; or (2) less pervious om the of input parameters 


in uniform profiles the capillary drive term in the Green and Ampt ‘equation — 
was best estimated by the water entry value as suggested by Bouwer (4). In 
their study, the capillary drive term calculated from soil-water properties was 
consistently greater than that obtained from direct measurements. Aggelides — E 
and Young (1) claim that these errors stem from the approximations inherent — - 
4 in the derivation of the Green and Ampt equation from the moisture flow equation. a 
ea Detailed descriptions and analysis of the currently available methods of calculating 
; the capillary drive terms are beyond the scope of this paper. The reader should 
consult the references previously mentioned for further information, = 
_ Estimation of the input parameters for Case 2, in which the surface i 
is less pervious than the subsoil, is 1 not as $ straight forward as § for Case i 


a a negative pressure © develops in . the zone below the ‘surface layer and remains 
constant for a considerable depth. Tagaki’s findings imply that if the subsoil 
is assumed homogeneous, then the moisture content and hydraulic conductivity 
inthis zone must also be constant and less than the values at saturation. — 
_ anegative pressure would not develop in this zone. These assumptions of constant 
_ moisture content and hydraulic conductivity are inherent in the development 
of both the GAML and modified GAML equations. The equations, then — 
are directly applicable to Case 2, of which the problem of surface sealing is — 
probably the most important example. For Case 2, He and S,,, can be evaluated — 
using the aforementioned techniques for Case Variables K, and can 
be evaluated by considering steady-state infiltration. 
s For steady-state infiltration gravity forces dominate and the capillary forces” 
can be neglected. The hydraulic gradieni approaches unity and the infiltration _ 
rate, given by Eq. lla, approaches the hydraulic thse cnet in the wetted — 
a For steady-state infiltration, the flux through the restricting ‘surface im. ; 
‘must equal the flux. the subsoil transmission zone, 


in in which 0 = the hydraulic | head. The hydraulic gradient in the ae transmission 
zone under -state conditions tends to so 


i 
Boor K, 


epth of ponding on the soil surface is negligible om bax. 


the capillary suction immediately below the ‘surface layer; and 


—— 


| With a sou of the K-S relationship, and the known value of K,/L,, 
_K, can be determined from Eq. 206 by trial and error. The moisture content | 
"behind the wetting front in the subsoil zone corresponds to this value of “4 
is determined from the S-6 relationship. From these values, and S_, 

Alternatively, the model parameters can be determined directly by employing 
_ functional relationships between | capillary suction, unsaturated hydraulic con- 

_ ductivity and moisture content, such as those proposed by Campbell (6). . These 

methods are only valid if a log-log plot of the moisture retention curve is 

a Straight line, and tend to break down close to saturation. Care should therefore 

be exercised in using these techniques to calculate K, and 

S Surface Sealing. —The modified GAML equations eee a simplified m means 


tm, 


q 
: 
1 


_ In such ¢ cases, £ ' will be ‘essentially constant, and K, can n be determined using. 7 
a _ the otecementiones technique. Eqs. 15 and 17 can then be solved in on same 


100 


NUMERICAL SOLUTION 


sw) 


INFILTRATION RATE (cm 


FINE /COARSE —» 


FIG. 3.—Predicted and Infiltration Rates Versus Time for 
Sandy Loam) Over-Fine (Millville Silt Loam) and Fine Over-Coarse 


(L.=Sem 


@ 
‘Ifth 
in whic 
L, = th 
— 
| 


manner as the original GAML equations u using g the suggested Mein 
and Larson (16,17). That is, Eq. 15 can be solved directly for t, or by iteration 


In many cases, K, is ‘not constant. For example, Moore et al. (19) suggest 


in which K, = the conductivity of a stable, 
— lished surface seal; K, is the initial hydraulic conductivity; and a = a constant — 
- controlling the rate ‘of decay of the initial hydraulic conductivity. Incssion 
of Eq. 21, or any other relationship to describe the change in K, with ry 
produces complicated equations that are difficult to solve (K,, K,, and the _ 
moisture content in the subsoil transmission zone all change with time). Examining : 
these effects is beyond th the of this paper. 
ANavysis 


n Djavan (7) measured infiltration rates into two-layered soil columns consisting _ 


of Lewiston fine sandy loam (coarse) overlying Millville silt loam (fine). He 
also measured the K-S and S-6 relationships for both soil types. These data 
permit the modified GAML equations to be evaluated against numerical solutions 
of Richards’ equation and measured infiltration rates. The numerical technique 
used here is ; essentially the same as that described by Smith and Woolhiser 
(23). Two cases are examined; (1) Coarse over-fine; and (2) fine over-coarse — 
The infiltration rates predicted by the numerical method and the modified 
_ GAML equations are shown for both Cases | and 2 in Fig. a together with — 

- observed infiltration rates (7) for Case 1. The application rate was 8.0 onfh, 
gy the surface layer was 5 cm thick. The calculated input parameters for 
the modified GAML equations are shown in Table], 
For the coarse over-fine stratification the times to surface ponding agree 
“very well. Predicted infiltration rates for the numerical method and the modified 
 GAML model also agree ‘Teasonably \ well. The difference between the two <i 


TABLE 1.—Input Parameters. for Modified GAML Equations for Fine" and Coarse 


7 


Conductivity ‘ | centimeters per hour 


q 
a 
| 
| 
| 


due to the failure of the wetting front to conform to the | pleten -flow assumption 
inherent in the derivation of the GAML equations. Observed infiltration rates 
are higher than predicted at intermediate times (14 min-30 min), but are in 
very close agreement with the predicted at larger times (>30 min). a 
Pett of both the times to surface ponding and infiltration rates for 
the fine over- “coarse Stratification show close agreement, except between 12 
min and 19 min from the beginning of the event. The reason for this arises, 
‘agile. from the failure of the actual wetting front to conform to the piston-flow. 
4 model. In the modified GAML equations the piston-type wetting front encounters 
the interface between the two strata at 18 min, at which time the infiltration | 
rate drops to an approximately constant rate. The leading edge of the ; wetting 
as predicted by the numerical method, encounters the interface at 11 
‘Reesiens and so produces a gradual reduction in the infiltration rate to a more 
or less constant value after it passes through the interface. The most significant _ 
q _ result here is the fact that both the numerical method and the modified GAML — on 
equations predict this accelerated reduction of infiltration rates when the wetting 7 a 
front encounters the interface, though the latter predicts a ‘much m more tent 
reduction. In addition, the constant infiltration rates after the worn front 
TABLE 2.—input Parameters for GAML Equation for Ida Silt Loam Soil (16) 


Parameter 
Wed 
Total porosity = centimeters per 
Capillary drive at wetting front 7. "centimeters of water 


Saturated hydraulic conductivity K. centimeters perhour 


q 


- moves through the interface are almost identical. This tends to validate the 
_ method of calculating K,, S,,, and 46,, for fine over-coarse strata described 
earlier, and shows that the input parameters can be measured, and thus do 
Asseed and Swartzendruber (3) have experimentally evaluated Eq. 105 for 
two-layered soil profiles consisting of coarse over-fine stratifications. ‘They _ 
- evaluated the parameters in the equation by fitting to infiltration data measured 
in both uniform and stratified profiles. For coarse over-fine stratifications their 
results were good. For fine over-coarse stratifications, their results were generally — 
reasonable for predictions based on the parameter set obtained by fitting to 
stratified profiles. Predictions based on the parameter set “obtained by fitting ay 
to uniform profiles were poor at larger times. This is expected because wae 


parameters do not take into account the restricting effect of the fine surface 
layer, which prevents the coarser subsoil from reaching saturation, = 8 — 
= In order to illustrate the magnitude of the surface effects predicted by the — 
modified GAML equations, the following hypothetical situation is presented. 
The parameter set derived by Mein and Larson (16) from measured K-S and 
'S-6 relationships, for an Ida silt loam soil will be used. The input parameter 


® Ind In the ——— the conductivity of the surface and subsurface layers of soil 


| 
| 
| 


pow ta oail was assumed). The surface seal was ‘assumed are be 5 mm thick 

_ (8) and the initial moisture deficit and capillary drive at the wetting front were 
assumed to be the same in the surface and subsurface zones. Analysis was 
performed for the various combinations of situations resulting from two different 
initial moisture contents (6,, three different rainfall intensity ratios (J/K,), and ; 
four different conductivity ratios (K, /K vs in \ which ‘K,= = the saturated anal a 
in the subsoil zone. | 


Initial 


is 


“Calculated by Richards’ equation (16), = 


TABLE 4.—Time to Surface t,, in minutes 


~ 


=" 

Initial moisture ‘Conductivity <T 

contents” fatio, (K,/K, = 


(3) 
276.4 


202.7 
138.2 


7 7 Table 3 compares F, predicted by the GAML model (equivalent to the modified - 


7 -GAML model when K ;/K, equals one) to values calculated from Richards’ 
equation by Mein and Larson (16). The error, expressed as a percentage, ranges 
- from —33%-38%. Mein and Larson (16) compared predictions of the GAML — 
bed model to those of Richards’ equation for four different soils. They ol 
that the Ida silt loam yielded the worst results, classifying them only as ‘‘fair.”’ 
‘Table 4 shows the time surface ponding, 1 “predicted by the modified 
_GAML model. Time ¢, decreases with increasing rainfall i intensity ratio, lester 


conductivity ratio, and increasing initial moisture content. 


‘ q 
TABLE 3.—Predicted and Computed Cumulative Infiltration Volume at Time to Surface 
— Ponding for K,/K, = 
0.28 
| Intensity Ratio, 
= 
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tie infiltration as functions of time for WK. to 
+ é are shown in Table 5. The cumulative infiltration volumes for 6, equal | to ~ 
0.40 are also shown in Fig. 4, together with the cumulative infiltration volume ea. 
calculated from by Mein and Larson (16) equal 


TABLE 5. Cumulative Infiltration in centimeters, Versus Time for 11K, = = =4 


0.730 


— — richaros’ EQUATION 


TIME (hours) 


OF Figs. 5 and 6 are different ways of comparing ‘the computed infiltration cates : 
as functions of time. Fig. 5 clearly shows the effect that the conductivity ratio 
has on both the shape and position of the infiltration rate curve, as well as 
the time to surface ponding. In a similar \ way, Fig. 5 shows the effect of rainfall _ | 


0.388 | 0.648 0.862 | 
| 


ta) 


hr 
(hours) 
FIG. §.—Infiltration Rate Time for Four Ditforent Ratios (ida 


i= 


. 6.—Infiltration Rate Versus for Ratios (Id 


conductivity of a thin surface layer decreases the infiltrability significantly and 
distorts the shape of the infiltration rate curve at high conductivity ratios. The 
- cumulative infiltration volume is reduced by: 45% or more after 5 h for conductivity 


can have on “infiltration the hydraulic 


The GAML infiltration model was modified to include surface effects such _ 
- surface sealing. For infiltration of a constant intensity rainfall greater than 
a the saturated hydraulic conductivity of the soil (J > K,) into a soil of uniform | 
moisture content, the the governing can be in the following 


4 
q 
| 
| 


form: & F, in which F, is given by either of the 


ae Ls 


+ 


| 


\E+F- 8) 
+ 


of of ‘the modified GAML for case of a 
«constant surface layer hydraulic conductivity are presented. These 
ls illustrate the importance of considering surface effects in infiltration studies 
_ and show the of these effects. 
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capillary suction at ‘wetting front; 
equivalent time to infiltrate volume, F,, under ponded surface 

vertical distance from soil surface; 
initial moisture deficit; and 
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f = steady state; sie 
= conditions at instant ponding 
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_ In the management of an irrigation project, it is important that the flow at 


canal bifurcations be measured accurately, for only in this way can the efficiency 
of irrigation water use be calculated. Accurate flow measurements are also 
important in drainage system management since they | are the / only means of 
assessing a reliable water balance upon which drainage criteria. can be based. ‘nd 
3 In many irrigation or drainage systems, however, the available energy headloss 
over a discharge measurement structure is limited. For accurate flow measurement | 
in such channels, the long-throated flume is the most promising structure (2,3) 
because it can be tailored to fit the stage-discharge curve of the channel in ; 
which it is to be placed, it is relatively easy and cheap to construct (2) and, 
‘for large flumes, the critical submergence ratio can be as high as 0.95 (1,2,4,5,7). 
. _ To obtain more data on this critical submergence ratio, as a function of the 
degree of downstream expansion, we ran extensive laboratory tests on two 
trapezoidal flumes. We also tested the influence of the slope of the entrance ai 


In a concrete- lined section, with bottom width, 5, = 0.50 m, and 

side slope of 1 to 1, a trapezoidal long-throated flume was placed (Fig. 1). 

‘The flume-throat was made of polyvinyl chloride | (PVC) and had a length, . 
= 0.60 m, width, b = 0.20 m, and a side ‘slope ratio of 1 to 1. The bottom © 
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$j” /|050(series 8) 
flow 0.81 (series C) 


downstream 
transition 


brackers refer 


replaced model 


O50=b, 
SECTION A- tor series SECTION OVER FLUME THROAT section & A- series C “wai 


FIG. 1 of Model, in meters, 1 m = 3.2808 


of the ‘throat wa was 0. 15 m the bottom of the approach ‘channel 


3 For individual series of tests, we fitted this throat Section with the following © 7 
entrance and downstream transitions: 

are A concrete entrance transition, the bottom of which converged | to 2. D. 


. transition of the side slopes started at the same distance upstream of the 
throat entrance (0.30 m) as did the bottom transition. The sides thus converged © - 
_ 1 to 1.5. Two downstream transitions were fitted to this flume, the angles of 
divergence being 1:0 (Test Al) and 1:6 (Test A5). We used these tests to test 
the influence of the short (1: 2) entrance transition on the e discharge | coefficient © 
2. As basic flume in the test series es B, the entrance : transition was was changed a a 


_ that the bottom slope became | to 3 (sides | to 2.25). This flume was tested _ 
_ with the downstream transitions shown in Table 1. During test B4 we encountered 

- some difficulties with the measurement of a stable tailwater level h >. To solve a 
- this problem, we first placed a pressure tap for } h, and for test BS we saseall 

r 3. For the C-series, the bottom of the approach and tailwater canals were 

- paised by 0.15 m so that p, = p, = 0. In this series the flat-bottomed flume | 

- thus needed no transition on the bottom. The angles of divergence a of the 

- side wall transition are the same as used in the B-series. ape 3 


ss ‘measured by a magnetic flow meter that was calibrated, before 
s tests started, against a gaged volumetric tank. This calibration was repeated | 7 
three times in between the tests. As a result, the error in the measured discharge _ 
_ is less than 1%. Heads were measured by point gages that can be read to 
the nearest 0. 1 mm. Each head given in in this paper is the average of of five readings. 


cenuiting ta aly ine 


in in which Q= in meters per second (Q in cubic seconds if 
United States are used i in the equation); C. discharge coeffi- 
cient (dimensionless); b = bottom width of flume throat, in meters; m = = side 

_ slope ratio of flume throat, horizontal to vertical; g = acceleration due to gravity, 
9.81 m/s’; H, = upstream energy head with respect to bottom of flume throat, 
in meters; andy. = critical depth = the flume throat, in meters. The last mentioned 

is a function of H,, b, , and m. For the tested flumes m = 1; y.-values 
a - The discharge coefficient C, has a predictable value provided that no separation — 
Pp 
. mr: of flow occurs at the boundary in the entrance transition. With a curved entrance | 
_ transition, e.g., cylindrical surfaces with sufficiently large radii, flow separation 
. 
z — will not occur. Construction costs, however, are lower with a plane surface 
entrance transition. . After a study by Wells and Gotaas (8), it became common 
TABLE 1.—Downstream Transitions Test Series B 


practice to use which 7 convergence each plane wall does 

“Rot. exceed 1: 3 Telative to > the flume center line. 


converged 1:2 and the side walls” 1:15. “The Cy -values ‘for these flumes are 
‘shown in Fig. 2. Also shown are the C,-values for two flumes, B, and B,, 
S which the bottom entrance convergence is 1:3 and side wall convergence 

q 1:2.25, and for one flume, C6, whose flat bottom plus side wall convergence 
is 25. For A, /L- ratios greater than 0. the Cy ~values of the A- “series are 


of the 6-flume. All C, -values, are well within the range of values 
7 derived from literature by the first writer (3). This range, which has a 95% 
confidence level, is shown in Fig-2, 

¥ On the basis of the aforementioned data, there would seem to be little objection 

‘in practice to using plane transitions, provided that the convergence of each 

7 wall does not exceed 1:2 relative to the flume center line. A typical flow pattern 

along the flume boundary is shown in Fig. 3. 

_ If the streamlines at the control section in the flume throat are straight and 
i and the flow is critical, the C,-value is independent of a patterns 

= downstream of the throat. ‘For ‘Tatio’s s of #A,/L | rea 
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TABLE 2.—Values of Ratio as Function ot H, tor Trapezcidel Control with 
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streamlines become curved, resulting in increasing C, 
for higher values of A, af L. In flat-bottomed flumes, (p, = =P2 — the discharging — 


in 


—Typical Flow Pattern Along Flume (A1) bo boundary, Q= 0.160 m 3/5 


transition seems to have no significant influence on Wal 
The Tandom distribution of C,-values for low values of A, / L (<0. 03) can 


_ curvature is thus less than in flumes with elevated throats. This difference J 
2.—C,-Values as Function of H,/L and Entrance Transition 
attributed to viscous effects, to the Changing velocity distribution in approach 


sf 
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FIG. 4.—C,-Values as Function of Downstream Transition and H,/L 
oneal. and flume throat, . ‘end to minor errors in the determination of =?) and — 
us , all of which have a relatively great influence on C,-values if H,/L is small. — 


the water level downstream of the flume atthe ff 
control section (in the throat) is supercritical and the discharge through the 
= Q, depends on H, only. This flow is referred to as free or modular. 
With rising downstream water level, however, flow at the control section becomes 
™ subcritical and Q depends ; also on the downstream energy head, H, - The modular 
4 1 limit is defined as the value of the submergence ratio, H,/ H,, at which the 
real discharge deviates by 1% from the discharge calculated by the head-discharge 
(Eq. 1). For of H, Eq. | be expanded to 
in which f = = drowned flow reduction factor, te less than unity. 
7 Fig. 5 shows a typical curve giving f-values as a function of the submergence “. 
- ratio. The higher the value of the modular limit, the more the curve will move 
- into the top right corner of this figure. In fact, we do not recommend submerged | 
flumes f for measuring discharge and and for the following reasons: 
. Errors in the determination of both H, and H, propagate into H,/H,. 
In the example of Fig. 5, an error of about 3% in the submergence ratio causes — 


We 

q | 

&g 
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REQUIRED toss 
“the I- value to vary from 0.82-0. 94, In Practice the the error in H,/H, often exceeds 


3%. If in addition the submergence ratio and the modular “limit are higher, 

the f- ~value i is so inaccurate that the true discharge can only be estimated. _ ae 
_ 2. If a flume is submerged (H,/H, exceeds the modular limit) while the: 
discharge is constant and H, rises, the upstream head must rise greatly to allow 7 
the discharge to pass. This may lead to overstopping of the upstream embankment . 
De termining and handling two heads is expensive. 
If a flume is to be designed in ind modular flow range, the modular -timit must — 
be known. For oF tested flumes, we determined this limit Sd nl by — 


a Q- H, of a very tom 
in = a 
ee dimension depends on u); and u = a dimensionless power deseution 4 
on the cross section of the flume throat. And secondly by calculating the head _ 


at the calculated deviates by 1% from the actually measured 


101) 


= log 1. 01 x H ) 
For each flume we set the measured discharge at 0.010 m 3) s, 0.030 m’/s, 
0. 090 m */s, and 0.160 m ‘/s, and raised the downstream water level stepwise. | 


— 
1 Submer- 
q 
| 


TABLE eee —Modular Limits and Seow! Losses for Tested Flumes as Function of 


6. 0253 | 0.0549 0.0963 


0.9002 
0.0154 


= 35. 3145 sec-ft, and im= = 3. 2808 ft. With flumes Bl, B2, and B3 


2. 43 m and for the other flames 3.83 m (see Fig. 1). 


| 


“a 
Ul arge, in cubic meters per é 
Tested ~ 0.010 | 0030 | 0090 | 0.160 
re h h 
AH 
— H,/H, 0.7812 
H, | 0.79890 | 0.8354 | 0.8138 | 0.7391 
h,/h, | 0.8176 | 0.8548 | 0.8136 | 0.7717 
0.8230 | 0.8607 | 0.8246 | 20.7901 
| 0.0214 | 0.0483 | 0.0745 
‘| 0814s 0.8567 0.8503 0.8216 
| aH | | 0.0385 | 0.0584 
| 08361 0.8757 | 0.8705 | 0.8416 
B6 H /H, 0.8414 0.8874 0.8796 0.8543 
| 07873 | 0.8232 | 09030 | osis 
AH | 0.0166 | 0.0228 | 0.0309 | 0.0533 
0.0152. 0.0231 | 0.0321 te 0.0501 
yh, | 0.8271 | 0.8702 0.8528 | 0.8654 
| 
| 
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and the matching h, . Table 3 lists the modular limits for each flume and pan a 


terms of both water levels, h,/h, and energy heads, H,/H, 


Whether flow through the flume is modular or ‘submerged | can | usually be 


_ from the way the jet leaves the flume throat. For modular flow the flow 
pattern is almost symmetrical with respect to the center line of the el 
If the flume is submerged, the jet moves either to the right or to the i 
as shown in Fig. 6. The side to which the jet moves can easily be influenced _ 
_by disturbing or changing the flow arent in the oe Via channel to the om 


With flumes in earthen channels, unstable jet can erode. the 


As can s can be seen seen from Table 3, the modular limit for the tested trapezoidal 
flume throat varies from about 0.70-0.90. Some general conclusions that can 
bedrawnfromthistableare; 

1. _ Downstream transitions with a ratio of expansion of 1:1, 1:2, and (to s some 

extent) 1:4, raise the modular limit only slightly. 

_ 2. Flat bottomed flumes, because of less streamline curvature, have a bicher 
— limit than flumes with elevated throats, 


; * _ 3. For the higher discharge range the modular limit is improved considerably _ 


for flumes B3-B6. For flat-bottom flumes (C3-C6) the modular limit varies 
only from 0.85-0.90. This is because most of the energy is lost in in the eddy 


irom 
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_ Although the aforementioned conclusions give some support in estimating 

_ the modular limit of a (trapezoidal) long-throated flume, we still needed a generally e 
_ applicable method of estimating the modular limit, and thus the energy head 

loss over a flume. In our search for such a method, we divided the loss of 


4 


(energy) head into: (1) The head loss between the upstream head measurement 
section and the control section in the flume throat; (2) the losses due to friction 
between the control section and the downstream head measurement section; - 


throat bottominm 
— 


-profile I 


Ad —The energy losses over the section upstream of were 
given by the first writer (2) whim od fea 


for the control section in which flow i is critical. In Eq. 8: a = the energy coefficient — 
_ that take takes into account nt the nonuniform distribution of velocities across ss the s section; 


and on the pressure distribution. The -value can be calculated by A 
Table 2, after which v, can be derived from the measured — mor wee - 


= 
| 
ue tO energy CONVErsion OVE! = 
i 


For values of H,/L< 0.4, the control section are are 
* ‘Straight and the velocity distribution is uniform, as is shown in Fig. 7. Thus, | 
it may be assumed that a = 1 and B = 1. Eq. 7 then gives a good estimate > 
7 of the energy head loss upstream of the control section, as is shown in Fig. 
If the ratio, H H,/L, increases, Streamlines at the control section become 
increasingly curved, thereby influencing the values of a, B, and Cx However, 4 
- within the limits of application of H,/L for long-throated flumes (2), 0.1 <_ 4 
A, /L < 1.0, Eq. 7 gives a good estimate of the energy head loss upstream — 
Ad (2).— —Although flow is nonuniform in th the downstream transition, we 
calculated the energy losses due to friction by applyi ing the equ equation yn of Manning 


to three reaches: (1) Reach of the fume throat downstream of the control 


= 


= 


Soe =——_- 


ei. 


FIG. 8. over Transition, (1 mm = 0.03937 in.) 
- 4 section, which approximates ( /3) L = 0.20 m; (2) length of reach of the actual | 
- transition of (bottom and) side walls; and (3) length of canal reach from transition 


- to measurement section of h, 2. (In our calculations, we used values of n = 
0.011 for PVC, n = 0.013 for plywood, and n = 0.015 for concrete. rae 


4 _ The total energy om due to friction, AH,, between the control section mand 


Table 4 lists the calculated values of AH,, together | with values of AH and — 
As can be seen in Table 4, the energy loss due to friction, H, — H, + AH,, © 
asa of the total energy loss over the flume, ‘AH, ‘increases 


— 
= 


with i increasing —" of the downstream transition. For a flat-bottomed flume 

with 1:10 transition (C6), this percentage is already larger than 35. Because 

of the relatively high flow velocities in the transition, the percentage of energy 
loss due to friction is larger in a flat-bottomed flume (C-series) than in a flume 
with an elevated throat (B-series). In the latter type, a gradual expansion can 
: Ad (3).—The losses due to energy con — on over the downstream transition 
. are assumed to equal the total energy losses minus all losses due to friction. _ 


TABLE 4. —Sasay Losses over Flumes | Flowing at their Modular Limit, in centimeters 


——— _ in cubic meters in cubic meters 
per row, in per secon Tested per secon 


4 
j 


Note: 1 m’/s = 35.3145 sec-ft, lom = 0.032808 ft. 


‘These losses were analyzed with the generally known n equation: _ ee a 
H, — SH, = 
in which € = the energy loss coefficient, a function - the pos of expansion 
of the downstream transition; v. = average velocity in that section of the throat — = 
= flow is critical; and v, = average velocity in the section in which h, 
We calculated values of & for each flume and plotted. ‘them in Fig. 9 as a 7 
function of the angle of expansion of the transition. As expected, the — 
— decrease with 1 decreasing —- of expansion. ‘They are, however, well abo 
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A 
| @ 0.290 | 047 | Cl | 0.24 0.42 7 
184 269 | cz | 173 254 
H, 027 043 | C2 | 025 | 043 
AH, B2 | 007 0.07 0.18 
AH 148 | 214 | C4 2.31 
A, 0.18 | 0.46 | 031 
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the values of (6) for transitions with subcritical flow both th upstream 


-flumes t to be analyzed and é-values to be calculated. We could, however, use 
_ data from Blau (1), Engel (4), Inglis (7), and Fane (5). The €-values we obtained — 
in this way are shown in Fig. 9. They conteapond sen reasonably with the €-values 


losses due to conversion of kinetic into potential 


_ energy < over ver the downstream transition, we recommend the use of the envelope 


of expansion 


serie C (p, =p, = 
flume KW 


ENGEL 10, 13, 17, 2214 
flume1 ® @ 
flume2 @ 

vd 
INGLIS 39, 60, 8314 
@ as 


angle of expansion a, in degrees 


FIG. 9. iil of — as Function of Angle of Divergence of Downstream Transition 


The total energy loss over the flume at the modular limit can now now be calculated 


y adding Eq. 7, 9, » and 10, 


= H, - — H,Cl/" + AH, + 


The modular limit limit can be found by dividing both ‘Sides Ss of this e equation | by 


g 
serie B (p, =p, 
4 
7 
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Be. 12 is a general « expression n for r the modular limit of any / long-throated flume 
and can also be used for the, hydraulically similar, broad- crested weir. arts = 
Procepure To Estimate Limit tobe 
To find modular of a flume discharging Q at an upstream 
h,, both sides of Eq. 12 have to be equalized by trial and error. The oui 
Determine the cross-sectional area of flow at the station where 
measured, and calculate the average velocity, v,. 
Calculate H, = h, + v 
3. Calculate the ratio, H,/L, and determine the C, -value. (see Fig. 2). 
4. Read the power u direct from the head- -discharge equation or from a 


plot of A, versus Q on double- €-log paper. 
6. Calculate y, at the control section as a function of H, and of the throat 
“an 7. Determine the cross sectional area of flow at the control section Lael ; 


the waterdepth, Yes and calculate the average velocity v.. = | 
Use Fig. 9 to find an é- ~value as a function of the ‘angle of 
2 (Gu)estimate the critical value of h, and derive 


equation to 1/3(L) of the throat, to the transition length, | and to the « canal ? 
up to the h, measurement section. section. 

13. Calculate H, = h, + 
(414, Calculate H,/H,. 

7 15. Substitute the values from items 5, 10, 12, and Ad into Eq. 12. 

If Eq. 12 does not match, repeat items 9- 
Once some experience has been acquired, Eq. 12 can be matched after two 

or three trials. Since the modular limit varies with upstream head, it is advisable _ 
to estimate the modular limit at both minimum and maximum anticipated heads 


and to check if sufficient head loss, H, — A, is available in both cases. 


a Long- throated flumes with a prismatic preted cross section are very a. 
for the accurate measurement of flows in irrigation and drainage canals. These 


: flumes can be tailored to fit almost any range of discharge while the head - 


loss they require over the flume to ensure modular flow is less” than that that of 


all other known structures 


_ Extensive laboratory tests on the influence of the upstream and ‘seiiiandiliin 

& transition on the value of the discharge coefficient have shown that upstream 
a plane transition can be used, provides the convergence of each wall does : 
= exceed 1:2 relative to the flume center r line. An analysis of modular limit © ~~ 


| 

| 

— | 
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1. Downstream transitions with a ratio of pte of 1:1, 1: :2, and 1:4 raise 
2. Flat-bottomed flumes, because of less streamline curvature, have a higher | 
limit than flumes with an elevated throat. ay 


the elevated throat, and because thee energy losses duet to friction i ina a flat- worl 
flume may approach 40% of the total energy loss, a downstream transition, 


; Downstream of a flat-bottomed flume an expansion only slightly improves the 


The equ equation (Eq. 12) and procedure to estimate the modular limit, as presented 
; in this paper, can be used for all long-throated flumes with a prismatic throat 
cross section, and also for the, hydraulically similar, broad-crested weir. — a 
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The following syn symbols are use are used i in 1 this paper: 


cross-sectional area of flow (L’); 


B water surface width (L); 7: 
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= be bottom width (L); 


f drowned flow reduction factor eels ' 
acceleration due to gravity (LT~*); 
total energy head with respect to throat level (L); somali =. 
piezometric head with respect to throat level (L); 
length of throat i in direction of flow (L); 


height of bottom hump (L); 
= power of head (dimensionless); 
= average flow velocity (LT™'); thy 
= water depth (L); 


angle of divergence in degrees; velocity distribution coefficient 


= total energy loss over flume (L); 
AH, = energy loss due to friction downstream of control (L); and = ae 


critical flow w conditions; 
1 measurement section; and 
downstream section. 
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peeer or on any paper presented at any Specialty Conference or other meeting, the Proceedings 
of which have been published by ASCE. Discussion of a paper/technical note is open to 
anyone who has significant comments or questions regarding the content of the paper/technical 
note. Discussions are accepted for a period of 4 months following the date of publication — 
of a paper / technical note and they should be sent to the Manager of Technical and Professional — 


be extended by a written request from a discusser. - 


7 Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period | ay 


‘The original and three copies of the Discussion should be submitted on 8-1 1/ 2-in. (220-mm) 
by 1l-in. (280-mm) white bond paper, typed double-spaced with wide margins. The length of 
a Discussion is restricted to two Journal pages (about four typewritten double-spaced 
of manuscript including figures and tables); the editors will delete matter extraneous to the 
subject under discussion. If a Discussion is over two pages long it will be returned for shortening. - 
All Discussions will be reviewed by the editors and the Division’s or Council’s Publications : 
Committees. In some cases, Discussions will be returned to discussers for rewriting, or they - 


may be encouraged to submit a paper or technical note rather than a Discussion. — ‘ined M, wy 


‘Standards for Discussions are the same as those for Proceedings Papers. A Discussion is — 


subject to rejection if it contains matter readily found elsewhere, advocates special interests, j 
is carelessly prepared, controverts established fact, is purely speculative, introduces personalities, te 
J 


or is foreign to the purposes of the Society. All Discussions should be written in the third 


person, and the discusser should use the term “‘the writer’’ when referring to himself. The 
author of the original paper/technical note is referred to as “‘the author.” 2 ste 
. Discussions have a specific format. The title of the original paper /technical note appears 
at the top of the first page with a superscript that corresponds to a footnote indicating the 
month, year, author(s), and number of the original paper/technical note. The discusser’s full 
name should be indicated below the title (see Discussions herein as an example) together with 
_ The discusser’s title, company affiliation, and business address should appear on the first 
page of the manuscript, along with the Proceedings paper number of the original paper/technical _ 
Rote, the date and name of the Journal in which it appeared, and the original author’s name. __ 
Note that the discusser’s identification footnote should follow consecutively from the original a 


paper/technical note. If the paper/technical note under discussion contained footnote numbers 
1 and 2, the first Discussion would begin with footnote 3, and subsequent Discussions would 
Figures supplied by the discusser should be por by senate, Starting with A. This also 
applies separately to tables and references. In referring to a figure, table, or reference that 
in the original paper / technical note use the same number used in the original. 
It is suggested that potential discussers request a copy of the ASCE Authors’ Guide to 
Publications o ASCE f more detailed information Preparation aud submission of 
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STORAGE Coerrictents FROM Grounp-Warer Maps* 
Closure by Brij M. Sahni® and Harbhajan 


_ The discussers are in - general “agreement with the writers with r regard to nn 
practical utility and potential of the method which uses readily obtainable data 

without incurring extra time and cost in drilling observation wells for preliminary _ 
assessment of aquifer parameters. The basic principle in writing the flow equation | 7 
in a simplified form in Refs. 4 and 8 is essentially the same as used by the @ 
writers, namely, | reducing the differential equation to a simple finite difference 
algebraic equation. However, it was the intent “of the writers: 
_ the assumptions involved first in writing a general flow equation Jade taedd 


in which S, = specific storage coefficient, as 
‘Hols 
and then in reducing it to “ee ae finally, in arriving at the finite difference, 
_ Itis is important to appreciate te the restrictions imposed bye each of the: mathematical 
= involved for, as discussed at length by the writers, the practical utility 
and validity of the method is dependent to a great extent on whether or not 
: all these restrictions are satisfied within acceptable limits in a given field situation. 
instance, on page 161, the discusser, Soliman, states ou 
sa Flow through aquifers may be considered a two-dimensional problem when — 
a: aquifer thickness is nearly uniform. For an aquifer strip of finite 
_ width, along a flow path, the piezometric head fluctuation, A, can be — 
‘considered as a function of distance x and time t — 
It should be noted here that even under these conditions stated by i: 
the flow may not be unidirectional, for instance, as also discussed in the original _ 


x paper, at points close to any natural or artificial sources and sinks. ih veel 


“June, 1979, by Brij M. Sahni and Harbhajan S. Seth (Proc. Paper 14658). 
: > s Postdoctoral Fellow, Dept. of Irrigation Water Management, ne Rice Research 
Inst., P.O. Box 933, Manila, 
*Prof. and Head, Dept. of Civ. Engrg. _G. B. Pant Univ. of Agr. Tech., 
India. 
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It is also important to note that in Eq. 7 of Soliman’s discussion, s and 
T cannot and should not be always taken as constants. For example, even 
if the aquifer thickness were uniform, the lithology of the aquifer material may — 
i change, or the aquifer nature may change from a nonleaky artesian to a leaky _ 
_ one within a certain distance. In the latter om the flow in the aquifer is a 
by vertical leakage through other strata. Under these conditions neither s and ~ 
_T may be constants nor the assumption of unidirectional flow - of Eq. Ay re) 
3 Soliman would be valid. Such situations in the field are not un uncommon. — ~-. 
_ Development of Eq. 4 by the writers clearly brings out the limitations ea 
on the choice of finite space and time elements As and At. A judicious choice 
of these parameters based on actual observations rather than on guessing is 
very essential for the success of determination of S using groundwater maps. 4 
‘Soliman (4) recommends a water level record from a single well and ground- water 7 
piezometric map drawn from a single set of observations. Though a continuous 
record of ground water level at a point would give more details regarding ~ 
fluctuations at that point, it cannot help verify whether As and At chosen — 
calculation of AJ/As and Ah/At do meet the requirements for validity of Eq. _ 
4or Eq. 8. The four piezometric contour ‘Maps drawn by the writers for the 


ine 


justified to assume Ah/At over an ‘arbitrarily chosen Ax to be the same 7 
that observed at one single point as was done in Soliman’s study. On the other — 
hand, following the procedure described in steps 2 and 3 by the writers (p. 
207) it is possible to make judicious selection of the finite difference values _ 4 
of s and ¢ to meet the > requirements of Eq. 4 as closely | as desired. 
Ideally, if the resources available ‘permit automatic recorder readings at a 
-number of points scattered over the study area, more detailed data could be i 
_ obtained making possible even preparation of daily piezometric contour maps. 
- However, where that is not possible and the observations have to be made tf 
manually at a large number of points spread over an extensive area, a weekly 


or semiweekly (or more frequent, depending on the manhours available) collection ik 


of piezometric data would be adequate. But one ground-water map prepared 

_ from one-time observations and thereafter observing fluctuations of ground water — 

— level only at a single point to obtain Ah/At for the entire area as obo 
by Soliman is certainly not appropriate. < 
Errata. —The following corrections should be made to the original 


Page 209, paragraph 2, "subheading: Should ‘read Case I (Fig. 3a) of 


Page 211, ‘paragraph 


Case (Fig. 4) 
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Crop RESPONSE INFORMATION FOR WATER INSTITUTIONS” 


1 Closure by D. Lynne® and Roy R. Carriker’ 


: The discussion by Hatchett, Stevens, and Vaux on the original paper expresses — 
concern that the writers misunderstand the economic concept of scarcity, that 
. the writers fail to appreciate difficulties associated with development and use ; 
of crop response information, and that the writers advocate emulation of markets 
by government agencies ‘‘when markets themselves can accomplish the job.”” 
However, these concerns are not substantiated, aaa by the int intent of t the writers co 
or by the content of the article. 
_ The need to distinguish ‘ “‘scarcity”” and ‘‘shortage,”’ as expressed by the | 
od discussers reflects primarily a semantic difficulty. ‘‘Shortage”’ need not be 
_ temporary, self-correcting, or relevant only in the context of markets. Nonmarket : 
mechanisms for resource allocation do exist, and shortages can | occur and persist — - 
in the context the nonmarket allocation of resources. In any "case, 
“ existing terms of trade. . .”’ need not c carry ry short- term or long- term connotations, — 


need not presume a particular form of organization for resource allocation. 


The relevant point is that a decision framework of some kind is needed whenever 
resources are limited, demands appear unlimited, and choice is necessary. Making 
an issue of the definition of scarcity, and the need to distinguish ‘‘shortage”’ 
. as a separate concept, invites protracted debate over a point which is not at 
-all crucial to the subject of the article or to the context within which the subject 
_ Estimation of water demand functions is costly and complex. The discussers 
‘eaughasize this point and provide a useful clarification of the manner in a 
Ta in Eq. 2, is derived. This Process for calculating the MVP of 1 water is 


invariant) and competitive market conditions. The long run marginal value of 

irrigation water, as depicted by in Eq. 1, is not simply P, (0 Y/aW), however. 
Rather, r,,, in Eq. 1 represents the marginal value of irrigation water as all 4 

ri other inputs and water are varied simultaneously along the long-run, optimal 
expansion path. The discussers’ description of = P,(@Y/ is ‘Strictly 

_ The term r,,,, is referred to advisedly as Sea . marginal ability to pay . 
and not as ‘ “marginal willingness to pay,”’ since the demand for irrigation water 
is a derived demand—derived from the demand for the products of isrigated 

_ agriculture. ‘‘Marginal willingness to pay”’ is more appropriate with ‘reference 4 
to measures of the utility-based demand of consumers. 
‘The discussers ; assert that crop “response information is difficult and costl 
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to obtain, and that water allocation ir institutions are unlikely to have the budget 
and d expertise to acquire such | data. Nine this extent the discussers apparently — 


2 net value for irrigation water could be compared with estimates of net value 
in other uses. If there were a highly competitive market for water, the ‘net 
institutional setting does affect water The effects of price ‘support 
~ rograms would be reflected in the demand models through higher levels = 
product price (p in Egs. | and 2). Restrictions on production practices ioted j 
* affect the level of, and nature of, the inputs used and the character of the 
underlying production functions. Restrictions on crop type simply eliminate some 
of the possible water demand functions. All of these refinements can be handled 
within the conceptual framework provided by the writers. 
It is difficult to develop an ‘‘. . . abiding faith in the ability of government 
agencies to emulate markets administratively.’’ Certainly, nothing in the article 
. should be construed as such. However, administrative water law is a reality 
io in Florida. The purpose of the article, as ste stated in| the introduction, was to 
outline a conceptual basis for integrating crop-water response information into 
7 the decision-making process for water allocation with particular reference = 
the Florida system. That system does not currently involve any form of 
‘. The Model Water Code (29), which served a as the modei for the Florida Water 


Ref 62). Alternative approaches to allocation v was the subject of two 
by Kiker and Lynne (58,59), and included discussions of market approaches — 
as well as other approaches to water allocation. A thorough discussion of the — 
subject must await another forum. Three points are worth noting, however. 
First, , evaluation of alternative institutional fc forms, be they market, nonmarket, 
or some other form, requires some concept of what institutions are supposed 
to achieve and some criterion for evaluation. Proposals for market solutions, 
as well as other proposals, must be subjected to this criterion. "Analytical 
_ approaches to institutional design are the subject of several interesting articles 
4 [see, e.g., Ciriacy-Wantrup (56), Schmid (61), Randall (60), and Bromley (55)]. 
Second, the likelihood that a particular institutional alternative can be ¢ adopted 
depends upon how that alternative is assessed by those who are ‘empowered — 7 
‘0 effect institutional change. The socio-political climate may be such that markets i 
will never be allowed and used for the allocation of water, not withstanding 
the elegance of market approaches to resource allocation and the fact that > 
‘regulatory approaches entail a “| . costly and unnecessary burden’”’ to develop 
value information. ‘Finally, even if some form ¢ of market for water or water 
rights were designed and implemented, the importance 
= in each case to protect third parties and nonefficiency objectives -2e Wi 0 
_ still leave unresolved questions about how much the public sector niente a mn 
8 protect third party interests, or about the level at which such interests eet 
_be supported. These demand elements are real, even if they cannot be reflected 
in the market. Thus, the writers see a role for both market and no 7 
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Rut Densiry* 
5 er by George R. Foster* and Leonard J. see A. M. 


7 The authors contributed to erosion modeling for agricultural land by 
bf concepts from channel morphology and sediment transport theory for noncohesive 
_ channels. This discussion describes the | writers’ analysis of field measurements 
of rill erosion on agricultural lands. 
‘The authors assumed that agricultural soils are noncohesive when they suggested 
‘that critical shear stress for agricultural soils could be estimated from Shields" 
A: criterion (4). _ However, this estimate of critical shear stress for erosion of 
* agricultural soils is much too small, as the writers will show from analysis 
; 7 of field data (7). Critical shear stress was obtained (Fig. 2) from these data — 
by plotting observed erosion rate versus shear stress estimated from observed — 
data for discharge, velocity, slope, and unreported measured cross- -sectional 


R,= 0.504% ... 


in which R, = hydraulic radius, in feet; and A, = flow area, in square feet. a 
An estimate of critical shear stress from Fig. 2 is 0.06 
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Critical shear stress Shields « criterion is asa of 
D, s4 (Fig. 3), assuming a shear velocity of 0.23 fps [corresponding to the as 
; “level discharge of Meyer, et al. (7) study], and a specific gravity of 2.65 for 
Beye particles. Meyer, et al. (7) observed that aggregates greater than | 
- accounted for 14% of the sediment coming from the rills (7). From Fig. 
FF mm for D,, gives a critical shear stress of 0.013 psf, significantly less 
i than the observed value of 0.06 psf. The actual D,, of the soil was much smaller, 


which gives an even smaller critical shear stress from the Shields criterion. 
_ Therefore, the Shields criterion using D,, for primary particles making up oa 
agricultural soils like Russell silt loam significantly | underestimates critical shear 
3 stress for rill erosion. The Shields criterion has proven satisf actory for estimating | 
3 critical shear stress for transport of aggregates (specific gravity of approx 1.8) 
’ ont primary particles (specific gravity of 2.65) by flowinrills(5), 


soils: do not not armor. those | studied by N Meyer, et 


& SHEAR STRESS T (Ibs/ft 
FIG. 2.—Erosion Rate in Rills on Tilled Soil, Russell Silt L Loam, as | Function > 
al. did not during their tests. The shear stress required to detach 
hs particles i is greater then the shear stress required to transport coarse sand that 
productive midwestern ‘soils to armor during an annual cropping cycle. Land 
“Tre cultivated agriculture is usually tilled for seedbed preparation, -_ sometimes a 
of coarse e particles. If armoring g does « occur, more than year is usually needed 
a complete armor to develop. 


at rill ecosien. Poster (5) measured rill erosion on a cropland soil that had 
consolidated and had not been tilled for a year. The soil was tilled and the 
= measurements repeated. Rill erosion following tillage was f five times the rate 
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— the critical shear stress of the channel boundary. If this type. of oe 
_ is reached, rill erosion ceases. Rill erosion rates were essentially constant we, 
the duration of both Meyer et al. G, 

_ If any restricting layer is deep, the apparent equilibrium is one where the 
cross sections of rills erode to an equilibrium shape that erodes downward 
at an equilibrium rate for : ‘steady flow. If a rill reaches a nonerodible layer, 

it widens and eventually reaches a final width when rill erosion ceases. ll 


time required to reach a final width depends on depth to the nonerodible 


discharge, critical shear stress, and slope (6). 


& The authors’ channel geometric properties a, T,, and A, are functions of 
‘slope and soil (8). These functions could and should 
0 
= 


be included in Egs. 5, 6, 9, 11, 16, and 20-22. The authors arrived at a solution “A 
to their equations by assuming an equilibrium where shear stress of flow in 


_ the channel equals critical shear stress. This assumption is questionable, and 


‘. another concept is needed. Given channel morphology relationships for erosion —- 


; a cohesive soils, rill widths could be determined if Q, could be estimated — 

& ‘Discharge, Q.. in each rill is often an independent hydraulic variable and — 
is a function of rilling patterns which may be primarily deterministic rather — 


3 than random. Microtopography and macrotopography as affected by tillage, 
land form, and crop, influence rill frequency and discharge in each rill. a 
isc research on rill density should include studies on rill frequency and 


— FIG. 3.—Critical Shear Stress from Shields Criterion (4) for Shear Velocity of 0.23 sid : 
$q ‘ft/sec, Kinematic Viscosity of 1.05 x sq ft/sec, and Specific Gravity of 
a 
| 


MARCH 
“the fraction of ‘total discharge in each rill. “Also, study of rill morphology 
E ‘they ae as a function of properties of tilled, cohesive agricultural soils 
as they are affected by tillage, management, and cropping is another important 
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